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Abstract 

This paper has two parts: one about observational constraints, and the 
other about chemical evolution models. In the first part, the empirical 
differential metallicity distribution (EDMD) is deduced from three dif- 
ferent samples involving (i) local thick disk stars derived from Gliese 
and scaled in situ samples within the range, —1.20 <[Fe/H]< —0.20 
[Wyse, R.F.G., Gilmore, G., 1995, AJ 110, 2771]; (ii) 46 likely metal- 
weak thick disk stars within the range, —2.20 <[Fe/H]< —1.00 [Chiba, 
M., Beers, T.C., 2000, AJ 119, 2843]; (iii) 287 chemically selected G 
dwarfs within 25 pc from the Sun, with the corrections performed in 
order to take into account the stellar scale height [Rocha-Pinto, H.J., 
Maciel, W.J., 1996, MNRAS 279, 447]; in addition to previous results 
[Caimmi, R., 2001b, AN 322, 241; Caimmi, R., 2007, NewA 12, 289] 
related to (iv) 372 solar neighbourhood halo subdwarfs [Ryan,S.G., 
Norris, J.E., 1991, AJ 101, 1865]; and (v) 268 K-giant bulge stars 
[Sadler, E.M., Rich, R.M., Terndrup, D.M., 1996, AJ 112, 171]. The 
metal-poor and metal-rich EDMD related to the thick disk shows sim- 
ilarities with their halo and bulge counterparts, respectively. Then 
the thick disk is conceived as made of two distinct regions: the halo- 
like and the bulge-like thick disk, and the related EDMD is deduced. 
Under the assumption that each distribution is typical for the corre- 
sponding subsystem, the EDMD of the thick disk, the thick -|- thin 
disk, and the Galaxy, is determined by weighting the mass. In the 
second part, models of chemical evolution for the halo-like thick disk. 
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the bulge-like thick disk, and the thin disk, are computed assuming 
the instantaneous recycling approximation. The EDMD data are fit- 
ted, to an acceptable extent, by simple models of chemical evolution 
implying both homogeneous and inhomogeneous star formation, pro- 
vided that star formation is inhibited during thick disk evolution, with 
respect to the thin disk. The initial mass function (IMF) is assumed 
to be a universal power law, which implies the same value of the true 
yield in different subsystems. The theoretical differential metallicity 
distribution (TDMD) is first determined for the halo-like thich disk, 
the bulge-like thick disk, and the thin disk separately, and then for 
the Galaxy by weighting the mass. An indicative comparison is per- 
formed between the EDMD deduced for the disk both in presence 
and in absence of [0/Fe] plateau, and its counterpart computed for 
(vi) N = 523 nearby stars within the range, —1.5 <[Fe/H]< 0.5, for 
which the oxygen abundance has been determined both in presence 
and in absence of the local thermodynamical equilibrium (LTE) ap- 
proximation [Ramirez, I., Allende Prieto, C, Lambert, D.L., 2007, 
A&A 465, 271]. Both distributions are found to exhibit a similar 
trend, although systematic differences exist. In addition, the related 
empirical age-metallicity relation (EAMR) cannot be fitted by the the- 
oretical age-metallicity relation (TAMR) predicted by the model, and 
the reasons for this discrepancy are explained. 
keywords - galaxies: evolution - stars: formation; evolution. 

1 Introduction 

Abundance trends and the metallicity distribution within any Galactic sub- 
system are vital records to the formation and the evolution of the Milky Way. 
In particular, it can be wondered if each component underwent a separate 
chemical evolution both in space and in time, or some fraction of gas was 
transferred from an assigned reservoir into another one, implying violation 
of mass conservation. In the latter alternative, an additional question is if, 
and to what extent, simple models of chemical evolution can be used for each 
subsystem. To this respect, useful indications may be provided taking into 
consideration the G-dwarf problem. 

In general, the occurrence of a G-dwarf problem is related to the obser- 
vation of too few metal deficient G dwarfs (or of any other selected spectral 
class) with respect to that which could be expected from the Simple model 
of chemical evolution (e.g., Searle and Sargent, 1972; Pagel and Patchett, 
1975; Haywood, 2001). It has been established in different regions of the 
Galaxy: the solar neighbourhood (van den Bergh, 1962; Schmidt, 1963) and. 
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to a lesser extent, the halo (e.g., Hartwick, 1976; Prantzos, 2003) and the 
bulge (e.g., Ferreras et al., 2003). In addition, a G-dwarf problem has been 
detected in both bulge-dominated and disk-dominated galaxies (Henry and 
Worthey, 1999), which is consistent with the idea that the G-dwarf problem 
is universal (Worthey et al., 1996). 

The deficit of metal-poor stars (with respect to the prediction of the 
Simple model) may be interpreted in different ways, such as changes in the 
initial mass function (Schmidt, 1963; Adams and Fatuzzo, 1996; Bromm, 
2004; Bromm and Larson, 2004; Larson, 2005), inflow of unprocessed (Lar- 
son, 1974; Lynden-Bell, 1975) or processed (Thacker et al., 2002) material 
from outside, evolution with inhomogeneous star formation (Searle, 1972; 
Malinie et al., 1993; Caimmi, 2000, 2001a, hereafter quoted together as COC0; 
Caimmi, 2001b, hereafter quoted as COl; Oey, 2003; Karlsson, 2005; Caimmi, 
2007, hereafter quoted as C07), instantaneous recychng approximation (Hay- 
wood, 2001), or overlooking the thick disk population (Haywood, 2006). For 
additional alternatives and further details refer to earlier attempts (Pagel 
and Patchett, 1975; Pagel, 1989). 

Inhomogeneous (i.e. implying inhomogeneous star formation) models 
of chemical evolution succeed in both providing a solution to the G-dwarf 
problem and reproducing substantial scatter exhibited by the empirical age- 
metallicity relation (EAMR), with regard to the solar neighbourhood (e.g., 
Malinie et al., 1993; COO), halo stars (COl), and bulge stars (C07). The 
current paper aims to investigate if inhomogeneous simple models of chem- 
ical evolution are also consistent with the metallicity distribution in (i) the 
Galactic disk, conceived as made of three distinct subsystems: the metal- 
poor thick disk, the metal-rich thick disk, and the thin disk, and (ii) the 
Galaxy, conceived as made of five distinct subsystems: the above mentioned 
three, the halo, and the bulge. In any case, the related metallicity distribu- 
tion is deduced (by weighting the mass) from the data belonging to selected 
samples (assumed to be representative of the whole subsystem). 

Although the Galaxy is dominated (by mass and star number) by the two 
major subsystems, the bulge and the thin disk, still the contribution of the 
remaining ones, the halo and the thick disk, may be relevant in particular 
metallicity ranges. In addition, evidence for a halo-bulge and thick disk-thin 
disk collapse (Wyse and Gilmore, 1992; Ibata and Gilmore, 1995) implies that 
both the proto-halo and the proto-thick disk were more massive than their 
present counterparts. The empirical distribution of specific angular momen- 

"'^For this reference, two points may carefully be kept in mind, namely (i) values of a few 
parameters must be corrected as explained in Caimmi (2001b), Sect. 3, second paragraph, 
and (ii) the majority of figures do not correspond to the caption, as explained in the 
erratum (Caimmi, 2001a). 
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turn on the above mentioned subsystems, favours the idea of two different 
kinds of contraction (Wyse and Gilmore, 1992; Ibata and Gilmore, 1995). 
On the other hand, some evidence seems to exist on a continuous transition 
from an extended {R ~ 20kpc), pressure-supported halo, to an inner, flat- 
tened {R ~ 15kpc), rotation-supported halo (Chiba and Beers, 2000). The 
chemical abundance of thick disk stars suggests a similar history to those of 
metal-rich ([Fe/H]?a —1.3) halo stars, and the thick disk abundance patterns 
show excellent agreement with the chemical abundances observed in metal- 
poor bulge stars, suggesting the two populations were formed from the same 
reservoir at a common epoch (Prochaska ct al., 2000). 

The following samples have been taken as representative. For the thin 
disk, N — 287 chemically selected G dwarfs within 25 pc from the Sun, 
with the corrections performed in order to take into account the stellar scale 
heights (Rocha-Pinto and Maciel, 1996), hereafter quoted as the RM96 sam- 
ple. For the halo, = 372 kinematically selected subdwarfs (Ryan and 
Norris, 1991), hereafter quoted as the RN91 sample. For the bulge, N = 268 
K-giants in Baade's window (Sadler et al., 1996), hereafter quoted as the 
SA96 sample. For the thick disk, a fictitious sample has been made from the 
following sources: (i) a normalized metallicity distribution (total number of 
stars not reported therein) related to a volume complete sample of local thick- 
disk stars, derived from Gliese and scaled in situ samples within the range, 
-1.20 <[Fe/H]< -0.20 (Wyse and Gilmore, 1995), hereafter quoted as the 
WG95 sample, and (ii) N — Ad hkely metal-weak thick disk stars within the 
range, —2.20 <[Fe/H]< —1.00 (Chiba and Beers, 2000), hereafter quoted as 
the CBOO sample. Under the assumption that both WG95 and CBOO sam- 
ples represent to an acceptable extent the related thick disk population, and 
keeping in mind that they overlap within the range, —1.20 <[Fe/H]< —1.00, 
a fictitious sample oi N — 592 stars, hereafter quoted as the FS07 sample, 
can be derived for the thick disk. 

The oxygen abundance is deduced from iron abundance, using two al- 
ternative relations involving the presence (PP) or absence (AP) of [0/Fe] 
plateau for sufficiently low [Fe/H] values (COl, C07). Though more precise 
dependences have been found in recent investigations (e.g., Jonsell et al., 
2005; Fulbright et al., 2005; Garcia Perez et al., 2006; Melendez et al., 2006; 
Ramirez et al., 2007), related to different trends for different Galactic sub- 
systems, still they arc lying between the above mentioned alternatives, which 
shall be used to preserve comparison with earlier results (COl, C07). 

The value of solar oxygen abundance also affects the [0/Fe]-[Fe/H] rela- 
tion. Recent determinations point towards lower values (e.g., AUende Prieto 
et al., 2001; Sofia and Meyer, 2001; Asplund et al., 2004; Melendez, 2004), 
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but the question is still under debate (e.g., Landi et al., 2007; Socas-Navarro 
and Norton, 2007). A solar oxygen mass abundance, (^o)o = 0.0056, de- 
duced from AUende-Prieto et al. (2001), shall be used to preserve comparison 
with earlier results (COl, C07). 

Due to the above mentioned uncertainties, it has been preferred to deal 
with less recent EDMD determinations to preserve comparison with earlier 
work (COO; COl; C07). More recent, statistically more significant and less 
well scrutinised determinations e.g., the Geneva-Copenhagen survey (Nord- 
strom et al., 2004) for the local disk and the ongoing Hambourg survey 
performed by the Beers consortium (e.g., Prantzos, 2007) for the metal-poor 
halo, are related to iron abundance instead of oxygen abundance as in a se- 
lected sample of nearby stars (Ramirez et al., 2007), and for this reason shall 
not be used in the current attempt. 

The empirical, differential metallicity distribution (EDMD), if) = AN/ 
{NA(j)), log0 =[0/H], has been determined in earlier attempts for the RN91 
sample (COl) and the SA96 sample (C07), using both the PP and AP depen- 
dence of oxygen abundance on iron abundance. The same is done in Section [5] 
for the RM96 and the FS07 sample. In addition, a putative EDMD is derived 
for the thick + thin disk and the Galaxy, respectively, by weighting the mass 
of the related subsystems, in Section [3l A comparison with the predictions of 
simple models, involving homogeneous and inhomogeneous star formation, is 
given in Section H] and El respectively. The discussion and the conclusion are 
the subject of Section [H] and [3 respectively. 

2 The data 

2.1 The empirical differential metallicity distribution 

In dealing with simple models of chemical evolution, involving the assump- 
tion of instantaneous recycling, the predicted metal abundance has to be 
compared with the observed oxygen abundance (e.g., Pagel, 1989; COO; COl; 
C07). Unfortunately, oxygen is more difficult than iron to detect, and an 
empirical formula is needed, to express the former as a function of the latter. 

On this matter a clear dichotomy appears between authors who support 
a plateau in [0/Fe] for stars with [Fe/H]~ —1 (e.g., Carretta et al., 2000); 
and those who do not (e.g., Israelian et al., 2001a, 2001b). Further insight 
can be gained by referring to specific proceedings (edited by Barbuy et al., 
2001). 

Taking the above alternatives as limiting situations, the following formu- 
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lae (COl; C07) shall be used: 
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in presence or in absence of [Fe/H] plateau, respectively. The value used for 
the oxygen solar abundance, (.^0)0 = 0.0056=f0.0006, has been deduced from 
a number ratio, log^no/nn) = —3.31 =F 0.05, or Aq = log{no / uh) + 12 = 
8.69 =F 0.05, corresponding to the solar photosphere line at 6300 A(Allende- 
Prieto et al., 2001), assuming a hydrogen abundance, X = 0.71 (Anders and 
Grevesse, 1989). For further details refer to earlier work (COl). The assumed 
oxygen abundance agrees with recent results (e.g., Melendez, 2004; Asplund 
et al., 2004, 2005; Shchukina et al., 2005; Ramirez et al., 2007) in contrast 
with earlier determinations (e.g., Sauval et al., 1984; Grevesse et al., 1984; 
Anders and Grevesse, 1989) which yielded higher values {Aq ~ 8.9), but the 
question is still under debate (e.g., Landi et al., 2007; Socas-Navarro and 
Norton, 2007). 

A number of [0/H]-[Fe/H] dependences lying between those expressed by 
Eqs. ([I]) and ([2]), but different for different Galactic subsystems, have been 
derived from recent investigations (e.g., Bensby et al., 2004; Jonsell et al., 
2005; Fulbright et al, 2005; Garcia Perez et al, 2006; Melendez et al., 2006; 
Ramirez et al., 2007). To allow comparison with previous results related 
to halo subdwarfs (COl) and K-giants in Baade's window (C07), the solar 
oxygen abundance assumed therein (Allende Prieto et al., 2001) shall be used 
in the current attempt together with Eqs. ([T]) and ([2]), hereafter quoted as 
"in presence" and "in absence" of [0/Fe] plateau, respectively, with regard 
to sufficiently low metallicities, [Fe/H]~ —1. 

Of course, the validity of Eqs. ([T]) and ([2]) shall be restricted to an epoch 
where only pop. II stars were present, [Fe/H]< —4 say, when the formation 
of extremely metal-deficient ([Fe/H]< —4) but oxygen overabundant (e.g., 
Christlieb et al., 2002; Bromm and Loeb, 2003; Iwamoto et al., 2005; Frebel 
et al., 2006) stars had already been occurred in presence of pop. Ill stars. 

While observations are related to logarithmic number abundances, [A/II] = 
log(A/H) — log(A/H)0, models of chemical evolution deal with mass abun- 
dances, Za = M\/M, where A is a generic element heavier than He, Ma is 
the total mass under the form of element. A, and M is the total mass under 
the form of any element. The following relation (Pagel, 1989; Malinie et al.. 
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1993; Rocha-Pinto and Maciel, 1996; COO; COl; COT): 



log (j) = log 



(3) 



holds to a good extenll^. For a formal derivation refer to earlier work (C07, 
Appendix 1). 

Let Alog0 = A[0/H] = [0/H]+— [0/H]^ be a logarithmic, oxygen abun- 
dance bin deduced from A[Fe/H] by use of Eq. ([I]) or ([21). The related, oxygen 
abundance bin is: 

A0 = A+0 + A-0; A^0=|0-0^| ; (4a) 
-01 T + 

; = — ^ — ; (4b) 



= expio 



H 



where in general, exp^ defines the power of basis ^ and, in particular, exp 
defines the power of basis e, according to the standard notation. To take also 
the RN91 sample into consideration, bins in [Fe/H] equal to 0.2 dex shall be 
used (COl, C07). 

Following earlier attempts (Pagel, 1989; Malinie et al., 1993; Rocha-Pinto 
and Maciel, 1996), the comparison between model predictions and obser- 
vations shall be performed using the differential instead of the cumulative 
metallicity distribution, as it is a more sensitive test (Pagel, 1989) and al- 
lows direct comparison between different samples (Rocha-Pinto and Maciel, 
1996). The occurrence of a sensitivity error in performing observations, pre- 
cludes the use of a proper differential notation in dealing with the empirical 
differential metallicity distribution (hereafter referred to as EDMclfl). This 
being the reason why the bin length cannot be lower than the sensitivity 
error, and differential ratios, i.e. first derivatives, dN/ d0, must be replaced 
by increment ratios, AN/ A0, where A0 is the bin length, AN the number of 
sample objects within the selected bin, and is the total number of sample 
objects. 

Accordingly, the EDMD in a selected class of objects is defined as: 

^(0TA^0)=log^ ; (5) 

where the increment ratio, AN/A(j), used in earlier attempts (Pagel, 1989; 
Malinie et al., 1993) has been replaced by its normalized counterpart, AN/ 



^The oxygen mass abundance, Zq, has been defined as O in COO and COl, to simplify 
the notation therein. 

^It is quoted as EGD in COO and COl. 
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{NAcp), used in more recent investigations (Rocha-Pinto and Maciel, 1996; 
COO; COl; C07), to allow comparison between different samples. The uncer- 
tainty on AA^ has been evaluated from Poisson errors (e.g., Ryan and Norris, 
1991), as a^N = (ANY^"^, and the related uncertainty in the EDMD is (e.g., 
COl; C07): 



A^ifj = lip — tjj^l = log 



^ AN 



(6a) 



^ , ANTiANV/^ 
- NA4. • 

where ip~ — oo in the limit AN 1. For further details refer to earlier 
work (COl). 

The [Fe/H]-[0/H] dependence and corresponding mean fractional oxygen 
abundance, (p, and half bin width, A^(j), in presence of [0/Fe] plateau (PP), 
according to Eq. ([T]), and in absence of [0/Fe] plateau (AP), according to 
Eq. ([2]), respectively, are shown in Tab. [1] for the metallicity range of interest. 

As in earlier attempts (Pagel, 1989; Malinie et al., 1993; Rocha-Pinto and 
Maciel, 1996; COO; COl; C07), the bin sizes in normalized oxygen abundance, 
0, correspond to uniform bin sizes in [Fe/H], according to Eqs. ([3]) and (jl]) 
and Tab.[T], which implies non uniform bin sizes, A0 = 2A'''0, in the following 
tables and figures. 



2.2 Metallicity distribution in the thick disk 

Both the WG95 and CBOO samples are considered to infer the metallicity 
distribution in the thick disk. The total number of sample objects is not 
specified in the former the related abundance distribution is nor- 

malized to unity (Wyse and Gilmore, 1995). A fictitious value, = 39, is 
deduced from Tab. 3 therein, and shall be used in performing calculations. 
The metallicity range is -1.20 <[Fe/H]< -0.20 for the WG95 sample and 
-2.20 <[Fe/H]< -1.00 for the CBOO sample. The related EDMD, using 
Eqs. dSD-dn]), is listed in Table [2] both in presence and in absence of [0/Fe] 
plateau. 

The samples under consideration overlap within the range, —1.4 <[Fe/H]< 
— 1.0, which makes a renormalization be possible. The result is a fictitious 
sample: N = 592, and metallicity range, —2.2 <[Fe/H]< —0.2. The related 
EDMD is listed in Table and plotted in Fig.[T] (top panels) both in pres- 
ence and in absence of [0/Fe] plateau. The trend exhibited looks like its 
counterpart related to the Galactic spheroid (C07). Having in mind it is 
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Table 1: The [Fe/H]-[0/H] dependence and corresponding mean fractional 
oxygen abundance, 0, and half bin width, A^0, in presence (PP) and in 
absence (AP) of [0/Fe] plateau. 
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Table 2: The empirical, differential metallicity distribution (EDMD) in the 
thick disk, deduced from the CBOO sample {N — 46, top panel) and the 
WG95 sample (A^ = 39 deduced from a normalized distribution, bottom 
panel), both in presence (PP) and in absence (AP) of [0/Fe] plateau. 
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Table 3: The empirical, differential metallicity distribution (EDMD) in the 
thick disk, deduced from the FS07 sample {N — 592), both in presence (PP) 
and in absence (AP) of [0/Fc] plateau. Bottom and top panels correspond 
to halo-like and bulge-like thick disk, respectively. The related EDMD, cal- 
culated for the halo-like {N = 18) and the bulge-like {N = 574) thick disk 
subsample, is denoted as ips- 
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33 


0, 


.205 


-0.378 


-0.361 


0. 


,137 


-0.033 


-0.020 


0.102 


0.082 


23 


0, 


.129 


-1.062 


0.458 


0. 


,101 


-0.784 


0.745 


0.374 


0.198 


3 


0, 


.081 


-0.640 


0.887 


0. 


,074 


-0.429 


1.101 


0.257 


0.161 


5 


0, 


.051 


-0.440 


1.063 


0. 


,055 


-0.296 


1.240 


0.257 


0.161 


5 


0, 


.032 


-0.240 


1.298 


0. 


,040 


-0.162 


1.364 


0.257 


0.161 


5 
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Table 4: The empirical, differential metallicity distribution (EDMD) in the 
thin disk, deduced from the RM96 sample {N = 287 and N = 325.064 after 
corrections due to stellar scale heigth; Rocha-Pinto and Maciel, 1996) both 
in presence (PP) and in absence (AP) of [0/Fe] plateau. 



PP AP 



















1.422 


-1.412 


1.604 


-1.587 


0.295 


0.174 


4.110 


1.129 


-0.166 


1.180 


-0.309 


0.061 


0.054 


57.410 


0.897 


0.136 


0.868 


0.027 


0.048 


0.043 


91.490 


0.713 


0.225 


0.638 


0.149 


0.049 


0.044 


89.127 


0.566 


0.004 


0.470 


-0.038 


0.072 


0.062 


42.593 


0.450 


-0.064 


0.345 


-0.073 


0.089 


0.074 


28.951 


0.325 


-0.614 


0.226 


-0.576 


0.151 


0.112 


11.843 



something else a mere coincidence, the fictitious FS07 sample is conceived as 
made of two fictitious subsamples related to different subsystems: a halo-like 
thick disk (A^ = 18) and a bulge-like thick disk (A^ = 574), in the metallicity 
range, -2.2 <[Fe/H]< -1.4 and -1.4 <[Fe/H]< -0.2, respectively. The 
corresponding EDMD is listed as ijjs in Table [3] and plotted in Fig. [21 top 
and bottom panel, respectively, both in presence (left panels) and in absence 
(right panels) of [0/Fe] plateau. 

The results from recent attempts involving a more refined abundance 
derivation (e.g., Reddy et al., 2006; Ramirez et al., 2007) are not used in the 
current paper, in absence (to the knowledge of the author) of their counter- 
parts related to the halo, the bulge, and the metal-weak thick disk. 

2.3 Metallicity distribution in the thin disk 

As in an earlier attempt (COO), the RM96 sample (A^ = 287 and A^ = 325.064 
after corrections due to stellar scale heigth; Rocha-Pinto and Maciel, 1996) is 
considered in the metallicity range, —1.2 <[Fe/H]< 0.2. The determination 
of the EDMD is repeated here owing to a different value of the assumed solar 
oxygen and, in addition, the absence of [0/Fe] plateau is considered. 

The results, obtained by use of Eqs. are listed in Table H] both 

in presence and in absence of [0/Fe] plateau. The number of objects in 
each metallicity bin has been calculated as AA^ = ANq/ f + 6{AN)i, where 
AA^o is the sample number, / and S{AN)i are correction factors for the 
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stellar scale height and observational errors plus cosmic scatter, respectively 
(Rocha-Pinto and Maciel, 1996). Accordingly, the total number of objects 
reads = 325.064. The related plots are shown in Fig.[T] (bottom panels) 
both in presence (left panels) and in absence (right panels) of [0/Fe] plateau. 
The dotted vertical line marks the transition from halo to bulge/disk glob- 
ular cluster morphological type (Mackey and van den Bergh, 2005). For 
further details refer to earlier work (C07). The distribution could be bi- 
modal with the occurrence of two maxima, related to low and intermediate 
oxygen abundance, respectively. The former maximum, if real, takes place 
near the transition from halo to bulge/disk globular cluster morphological 
type, as for old halo globular clusters and bulge K-giants in Baade's window 
(C07). 

Metallicity distributions with similar features are shown by other samples 
of thin disk stars (e.g., Wyse and Gilmore, 1995; J0rgensen, 2000; Kotoneva 
et al., 2002; Reddy et al., 2003). It has been questioned that the maximum 
of the distribution could be lowered by selection effects (Haywood, 2001). 
In fact, using different criteria, the maximum rises up to about the solar 
abundance (e.g., Favata et al., 1997; Haywood, 2001, 2006), but no general 
consensus seems still to exist on the selection of sample objects. For this 
reason, the RM96 sample shall be considered in the current paper. 



3 Inferred metallicity distribution 

Let the procedure used for the Galactic spheroid (C07) be generalized to a 
system made of ns subsystems, where the EDMD related to a representative 
sample is available for each subsystem. The total mass is: 

ns ns /i^. 

M = T.M. = Y.^M ■ (7) 

i=l i=l 

where Mj is the mass of i-th subsystem. 

Let A^ be the total number of long-lived (i.e. life time longer than the age 
of the system) stars in the system, and AA^ the number of long-lived stars 
within a selected metallicity bin. The relative frequency, AN/N, reads: 

AAT _ ^ ANi _^Ni ANi 

where AA^j and Ni are the number of long-lived stars within a selected metal- 
licity bin and the total number, respectively, with regard to i-th subsystem. 
Then the relative frequency, AN/N, related to the system, is expressed as a 
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mean of the relative frequencies, ANi/Ni, related to the subsystems, weighted 
by the factors, Ni/N. 

Under the assumption of a universal initial mass function (IMF) for star 
generation, the number ratios, Ni/N, may safely be related to the mass 
ratios, Mi/M, as: 

Ni _ Mi 

'N~1I ' ^ ' 

for further details refer to earlier work (C07). 
The combination of Eqs. (E]) and ([9]) yields: 



g M,, AN 



1=1 



M N 



(10) 



where the relative frequencies, ANi/Ni, may be deduced from representative 
samples, which allows an evaluation of the relative frequency, AN/N. The 
related uncertainty may be obtained using the standard formula of linear 
propagation of error^ together with evaluation of Poisson errors (e.g., Ryan 
and Norris, 1991), ctan, = (ANiY/^. The result is: 



0"AN/N = 2^ 



i=l 



M N 



which may explicitly be calculated for assigned samples and mass ratios. 

The EDMD related to the system results from the combination of Eqs. (( 
(ED, and dlO]), as: 



, , ^ M, AN 
^ = log 



M NiA(j) ' 



log 



i=l 



^ AN/N 
'M,AN,T{AN,y/^ 
M NiA(j) 



(12a) 
(12b) 
(12c) 



which may be particularized to the cases of interest. 



3.1 The thick disk 

The following assumptions have been made in Subsection 12 . II for deriving the 
halo-like and the bulge-like thick disk EDMD: (i) the parent samples WG95 

^Linear propagation instead of quadratic propagation has been preferred to maximize 
the errors. The symbol, A(AA^)/A^, would be more germane to this respect, on the left- 
hand side of Eq. (jlip , but ctan /n is used to avoid confusion. 
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and CBOO represent to a comparable, acceptable extent the thick disk within 
the metallicity ranges, -1.4 <[Fe/H] < -0.2 and -2.2 <[Fe/H]< -1.0, 
respectively, to allow renormalization with respect to the common range, 
—1.4 <[Fe/H]< —1.0; (ii) objects belonging to the resulting FS07 sample, 
whose EDMD cannot be fitted by simple models of chemical evolution, come 
from different subsystems, whose EDMD can be fitted by simple models of 
chemical evolution, namely the halo-like ([Fe/H] < —1.4) and the bulge-like 
([Fe/H] > -1.4) thick disk. In this view, the EDMD related to the thick 
disk depends, via Eqs. (fT2l) . on the halo-like to bulge-like thick disk mass 
ratio, Mhk/^bk, which (to the knowledge of the author) is not available at 
present. 

It has been suggested that the local density of the metal-weak thick disk 
population ([Fe/H] < —1) represents less than one percent of that of the 
canonical thick disk (Martin and Morrison, 1998) and, in any case, it appears 
to be a minor constituent of the entire thick disk population (Beers et al., 
2002). It remains unclear whether the metal- weak thick disk population is 
properly considered an actual subsystem of the canonical thick disk with 
a distinct EDMD, or whether it corresponds, in reality, to the low-metal 
tail of the canonical thick disk EDMD (Beers et al., 2002). The former 
alternative has been chosen in the current attempt, according to the above 
considerations. 

Different values of the mass ratio, Mhk/^bk; make points related to 
halo-like and bulge-like thick disk (on the left and on the right, respectively, 
of the vertical dashed line in Fig.[Tl top panels) EDMD shift vertically one 
with respect to the other: increasing values make the halo-like side of the dis- 
tribution rise and the bulge-like side lower, and vice versa. The special case, 
Mhk/Mbk = Nnx/NBK = 18/574 ^ 0.03136, leaves the thick disk EDMD 
(Fig.dl top panels) unchanged, as plotted in Fig.[3](top panels, crosses). The 
cases, Muk/Mbk = 0.01 (triangles) and 0.10 (squares) are also represented 
for comparison, where the error bars are suppressed to avoid confusion. 

3.2 The disk 

The disk is conceived as made of three main subsystems, namely (i) the 
halo-like thick disk; (ii) the bulge-like thick disk; and (iii) the thin disk. 
Accordingly, the EDMD related to the disk depends, via Eqs. ( fT2i) . on the 
halo-like to bulge-like thick disk mass ratio, Mhk/A^bk, and the thick to 
thin disk mass ratio, Mkd/j^nd, both unknown at present. The effect of the 
former ratio on the EDMD has been discussed in Subsection 13. 1[ 

The cases, Mkd/M^-d = 0.10, 0.05, and 0.20, are represented in Fig. [3] 
(bottom panels) as crosses, triangles, and squares, respectively, where the 
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error bars in the last two alternatives are suppressed to avoid confusion. In 
any case, Mhk/^bk = 18/574 ^ 0.03136 has been adopted. 

3.3 The Galaxy 

The Galaxy is conceived as made of six main subsystems, namely: (i) globular 
clusters (GC); (ii) field halo stars (FH); (iii) field bulge stars (FB); (iv) field 
halo-like thick disk stars (HK); (v) field bulge-like thick disk stars (BK); 
and (vi) field thin disk stars (ND). As in C07, the following mass values are 
assumed: 

Mgc = 0.001Mio ; Mfh = O.lMio ; Mps = Mio ; (13a) 
Mkd = Mhk + Mbk ; Mfd = Mkd + Mnd = 5.8M10 ; (13b) 
Mg = Mac + Mfh + Mfb + Mfd = 6.901Mio ; Mio = lO^V© ; (13c) 

where the halo-like to bulge-like thick disk mass ratio, Mhk/Mbk, and the 
thick to thin disk mass ratio, Mkd/Mnd, remain undetermined and must be 
fixed. 

The case, Mhk/Mbk = 18/574 ^ 0.03136 and Mkd/Mnd = 0.1, is repre- 
sented in Fig.m (middle panels) and compared with its counterpart related to 
the Galactic spheroid (top panels, C07). Additional cases, Mhk/Mbk = 0.01 
(triangles), 0.1 (squares), with Mkd/Mnd = 0.1; and Mkd/Mnd = 0.05 
(diamonds), 0.2 (asterisks), with Mhk/Mbk = 18/574; are also presented 
in Fig.m (bottom panels). The related error bars are not shown to avoid 
confusion. 

It is apparent that the effect of changing the halo-like to bulge-like thick 
disk mass ratio, Mhk/Mbk, and the thick to thin disk mass ratio, Mkd/Mnd, 
within reasonable ranges, yields an uncertainty on the EDMD which is of the 
same order as the Poisson error. Accordingly, from this point on the choice: 

Mhk/Mbk = 18/574 ^ 0.03136 ; Mkd/Mnd = 0.1 ; (14) 

shall be assumed as reference case. 

The EDMD in the Galaxy looks like its counterpart related to the Galac- 
tic spheroid (halo + bulge), which is bimodal. It is apparent that (closed 
or open) simple models of chemical evolution cannot provide a satisfactory 
explanation to a bimodal EDMD, and a different model is needed (C07). In 
this view, let the main Galactic components i.e. the halo, the bulge, the thick 
disk, and the thin disk, be supposed to undergo distinct chemical evolutions. 
Let the thick disk be supposed to be substructured into two subsystems, the 
halo-like and the bulge-like thick disk which, in turn, undergo distinct chem- 
ical evolutions. The related theoretical differential metallicity distribution 
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(hereafter referred to as TDMD) is calculated in a selected spectral class 
of long-lived stars, and the resulting TDMD is compared with its empirical 
counterpart. To this aim, simple models implying both homogeneous and 
inhomogeneous star formation shall be used, as in earlier attempts (COl; 
C07). 



4 Homogeneous simple models 

Homogeneous simple models with (inhibiting star formation) gas (Hartwick, 
1976) have widely been discussed in earlier attempts (COl; C07), and the 
interested reader is addressed therein for further details. What is relevant 
for the current paper, shall be reported in the following. 

Homogeneous simple models of the kind considered imply four main as- 
sumptions, namely (i) instantaneous recycling; (ii) homogeneous star forma- 
tion; (iii) universal power-law IMF; and (iv) gas inhibition from forming stars 
at a rate proportional to the star formation rate. Accordingly, the TDMD is 
represented as a straight line (e.g., Pagel, 1989; COO; COl; C07): 

dA^ 

^(0) =log^^ = a0 + 6 ; (15) 



and the explicit expression of the coefficients, a and b, reads: 

1 (^0)0 . 

In 10 p" ' 



(16) 



fe = logf^^^)-a0.; (17) 

where {Zo)q is the solar oxygen mass abundance, p" the effective (oxygen) 
yield, /i the (allowing star formation) gas mass fraction, and the indices, o 
and /, denote the beginning and the end of evolution, respectively. 
The combination of Eqs. (fT6!) and (ITTl) yields, after some algebra: 

^ = 1 + InlO aexpio(-a0o - ft) ; (18) 

and the condition, fij/fio > 0, translates into: 

6> -a0o + log(-lnlO a) ; (19) 
which makes a restriction on (otherwise acceptable) linear fits to the EDMD. 
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The oxygen mass abundance, Zq, to a large extent, may be related to the 
gas mass fraction, n, as (e.g., Hartwick, 1976; Pagel, 1989; COO; COl; C07): 

Zo-{Zo)o = ~p"\n-^ ; (20) 

P"^^; (21) 

where p, is the true (oxygen) yield and k, is the ratio of (inhibiting star 
formation) gas mass fraction to long-lived star and stellar remnant mass 
fraction. 

It can be seen that a simple model with (inhibiting star formation) gas, 
lower mass limit of long-lived stars, rrimf, and inhibition parameter, k, is 
equivalent to a simple model with the same values of independent parameters 
except the above mentioned two, (m^/)! < frimf and Ki = < k; and vice 
versa. Accordingly, the related (true) yield, pi, has the same value as the 
effective yield, p", expressed by Eq. ( 1211) . For further details refer to earlier 
papers (COl; C07). 

For the models under consideration, the (inhibited from forming stars) 
gas mass fraction, D, is related to the long-lived star and stellar remnant 
mass fraction, s, as (e.g., Hartwick, 1976; COl; C07): 

D{t)-Do = K[s{t)-So\ ; (22) 

and mass conservation reads: 

/i(t) + s{t) + D{t) = 1 ; (23) 

where /i is the (allowing star formation) gas mass fraction. The combination 
of Eqs. (I22D and yields: 

— Tf^^ (24) 
D-D„^<!^; (25) 

finally, the related EDMD can be determined following a standard procedure 
(e.g., Pagel and Patchett, 1975; Caimmi, 1981), as performed in COl. 

With the aim of applying the above results to the chemical evolution of 
the Galaxy, the value of the true normalized yield, p/{Zo)q, the long-lived 
star and stellar remnant mass fraction in a star generation, a, and the lower 
mass limit of long-lived stars, rrimf, shall be determined from a linear fit to 
the EDMD in the disk solar neighbourhood (COO) with regard to a solar 
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oxygen mass abundance, (^0)0 = 0.0056 (COl), and inhibition parameter, 
K = 0. Positive or negative values of k, inferred from linear fits to the EDMD 
in a specified Galactic subsystem, would imply inhibited or enhanced star 
formation, respectively, in comparison to the disk solar neighbourhood. 

Two extreme values of the power-law IMF exponent, p, shall be consid- 
ered, namely: (i) p = 2.9, which is a fit to the IMF determined by Scalo 
(1986) or Miller and Scalo (1979), for m ^ (Wang and Silk, 1993), and 
provides a good approximation also in terms of oxygen production (Wang 
and Silk, 1993), and (ii) p = 2.35, which coincides with the IMF determined 
by Salpeter (1955). More recent fits to the Scalo (1986) IMF yield a gentler 
slope, p = 2.7 (e.g., Weidner and Kroupa, 2005). 

For fixed values of the normalized yield, p/Zq, and fractional oxygen 
abundance at the beginning and at the end of evolution, 0o and 0/, a change 
in the IMF slope leaves all the output parameters unaltered except the lower 
mass limit of long-lived stars, m^/, the mass fraction of a star generation 
locked up in long-lived stars and stellar remnants, a, and the efficiency of 
star formation rate, C, where the product, aC, is also left unchanged (COl, 
Appendix C, Theorem 2). Accordingly, values of the above mentioned pa- 
rameters, related to IMF exponents within the range, 2.35 < p < 2.9, lie 
between their counterparts related to p = 2.35 and p = 2.9. 

The characteristic stellar masses shall be taken as in earlier attempts 
(COO; COl; C07), namely: upper star mass limit, ruMf = 60; neutron star pro- 
genitor lower mass limit, m^^ = 9; long-lived star upper mass limit, rfimr = 1; 
oxygen synthesising star lower mass limit, rumro = 10.5; typical neutron star 
mass, m„s = 1.5; typical white dwarf mass, = 0.6; where m = m/m©. 
For references see e.g.: niMf, (Wang and Silk, 1993); m^s (Prantzos, 
1994); rutr (Prantzos and Silk, 1998); nimro (Pilyugin and Edmunds, 1996). 
On the other hand, the lower mass limit of long-lived stars, m^/, results 
from the model. 

Some parameters which remain fixed are listed in Tab. [5], where the in- 
dices, 2.9 and 2.35, denote the value of the power-law IMF exponent used 
in computing the corresponding quantities. The initial and final normalized 
oxygen abundance assumed for the halo-like thick disk, (pou and ipfu, and 
the bulge-like thick disk, and (pjB, are <^oh = 0.03; 0/h = 4>oB = 0.16; 
0/B = 0.80. The thin disk shall not be considered, as a large fraction of the 
related EDMD cannot be fitted by homogeneous models (see Fig.[T], bottom 
panels) . 

For an assigned subsystem, different models can be obtained by chang- 
ing a single remaining input parameter, chosen to be either the normalized 
effective yield, p"/{Zo)q, or the slope of the TDMD, a, which conforms to 
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Table 5: Values of parameters which remain fixed for different simple homo- 
geneous models where star formation is inhibited. Star masses are normal- 
ized to the solar value, m = m/m©. The indices, 2.9 and 2.35, denote the 
value of the power-law IMF exponent used in computing the corresponding 
quantities. 



p/{Zo)q ■ 10 


7.3722 


{mmf)2.9 ■ 10 


3.4235 


{mmf)2.35 ■ 10^ 


6.9136 


"2.9 ■ 10 


7.3666 


"2.35 ■ 10 


8.9104 


{Zo)e ■ 10^ 


5.6 




1 


So 





Do 






Eqs. f ll5p and ( fT6l) . A few cases are listed in Tab. [6], related to the halo-like 
thick disk (HK) and the bulge-like thick disk (BK), respectively, where 
represents the mean oxygen abundance (normalized to the solar value) of 
stars at the end of evolution. Model HK2 is very close to the limiting case, 
fif = 0, defined by Eqs. ( fTSj) and f fT9l) . Values of parameters related to inho- 
mogeneous simple models which fit the linear part of the EDMD in thin disk 
(DN), plotted in Fig.[I] (bottom panels), are listed for comparison. 

It is apparent that halo-like and bulge-like thick disk models demand (to 
a different extent) inhibited star formation to maintain (i) universal power- 
law IMF, and (ii) true normalized yield, p/{Zo)q, unchanged with respect 
to a value deduced from an acceptable fit to the linear region of the EDMD 
in the thin disk solar neighbourhood (COO). More precisely, about 73-80% 
and 49-71% of the initial halo-like and bulge-like thick disk gas, respectively, 
has to be inhibited from forming stars, for providing an acceptable fit to the 
related EDMD. It results in a normalized effective yield, p"/{Zo)q, about 
three-four times lower in the halo-like thick disk with respect to the bulge- 
like thick disk. It is worth remembering that the above difference raises to 
about one order of magnitude for the halo and the bulge, and enhanched 
(instead of inhibited) star formation is demanded for the latter. For further 
details refer to earlier work (CO 7). 

The TDMD deduced from models HK1-HK2 and BK1-BK2 is represented 
in Fig. [5] and compared to the corresponding EDMD in connection with halo- 
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Table 6: Values of parameters related to homogeneous simple models where 
star formation is inhibited, corresponding to linear fits to the empirical dif- 
ferential metallicity distribution (EDMD) in halo-like (HK) and bulge-like 
(BK) thick disk, plotted in Fig. [21 The mean oxygen abundance (normalized 
to the solar value) of stars at the end of evolution is denoted as (j). Model 
HK2 is very close to the limiting case, /x/ = 0, defined by Eqs. (fTSl) and (fTOll . 
Values of parameters related to inhomogeneous simple models which fit the 
linear part of the EDMD in thin disk (DN), plotted in Fig. [T] (bottom panels), 
are listed for comparison. The definition, Y = p"/{Zo)q, is to save space. 





HKl 




HK2 




BKl 




BK2 




DNl 




DN2 




Y 


8.3518 E- 


-2 


6.8285 E- 


-2 


3.6191 E- 


-1 


2.0508 E- 


-1 


4.4233 E- 


-1 


4.4233 E- 


-1 


—a 


5.2000 E- 


-0 


6.3600 E- 


-0 


1.2000 E- 


-0 


2.1176 E- 


-0 


9.8183 E- 


-1 


9.8183 E- 


4 


b 


1.3170 E- 


-0 


1.3885 E- 


-0 


8.1340 E- 


-1 


1.3311 E- 


-0 


1.0210 E- 


-0 


9.7643 E- 


-1 


K 


7.8271 E- 


-0 


9.7962 E- 


-0 


1.0370 E- 


-0 


2.5947 E- 


-0 


0.0000 E- 


-0 


0.0000 E- 


-0 




1.7355 E- 


-1 


7.1092 E- 


-2 


2.7260 E- 


-2 


1.7656 E- 


-2 


2.9047 E- 


-1 


1.2799 E- 


4 




9.3627 E- 


-2 


8.6040 E- 


-2 


4.7753 E- 


4 


2.7327 E- 


-1 


7.0953 E- 


-1 


8.7201 E- 


1 


Df 


7.3283 E- 


-1 


8.4287 E- 


-1 


4.9521 E- 


-1 


7.0907 E- 


-1 


0.0000 E- 


-0 


0.0000 E- 


-0 





8.2804 E- 


-2 


8.4469 E- 


-2 


6.2538 E- 


-1 


4.9020 E- 


-1 


5.9389 E- 


-1 


5.0003 E- 


-1 
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like thick disk (top panels), bulge-like thick disk (middle panels), and thin 
disk (bottom panels) stars, both in presence (left panels) and in absence 
(right panels) of [0/Fe] plateau. 

Though the linear fit to the thin disk EDMD is derived using inhomo- 
geneous models (COO), still it can be related to homogeneous models where 
a = -0.98183; ^"/(^o)© = 0.44233; and b = 1.0211 (PP), b = 0.97643 (AP). 
The dashed line in top panels corresponds to model HK2 in Table [HI which 
is very close to the limiting case, /i/ = 0, defined by Eqs. (ITSl) and ( fT9ll . 
and model HKl is represented by the less inclined (in absolute value) full 
line. Accordingly, the linear fits lying between the dashed line and the more 
inclined (in absolute value) full hne, must be excluded. 

For the halo-like thick disk, homogeneous simple models in the presence 
of inhibited star formation provide an acceptable fit (though in a restricted 
range) with values of input parameters as listed in Tab. [5] and 0oH = 0.03; 
0/H = 0.16; 0.068 < p"/{Zo)q < 0.084, both in presence and in absence of 
[0/Fe] plateau. 

For the bulge-like thick disk, homogeneous simple models in the presence 
of inhibited star formation provide an acceptable fit with values of input pa- 
rameters as listed in Tab. [5] and 0oB = 0.16; 0/b = 0.80; 0.205 < p" /{Zo)q < 
0.362, both in presence and in absence of [0/Fe] plateau. 

For the thin disk (represented for comparison), homogeneous simple mod- 
els provide no acceptable fit with values of input parameters as listed in 



Tab.Eland (poN = 0.325, 0/n = 1-58, (PP); (poN = 0.226, = 1-85, (AP); 
p7(Zo)0=p/(Zo)0 = 0.44233. 



Further inspection of Fig. [5] shows a marked deficiency in the number of 
stars observed below a threshold, [Fe/H]^ —1 for the bulge-like thick disk 
and [Fe/H]?« —0.8 for the thin disk, which is close to the transition from halo 
to bulge/disk globular cluster morphological type. In other words, a G-dwarf 
problem seems to exist for both the bulge-like thick disk and the thin disk. 

Let dNi be the predicted total number of long-lived stars within the i-th 
subsystem [1 < i < ns) belonging to a system, related to a fractional oxygen 
metal abundance, =F d0/2, with regard to a simple homogeneous model 
of chemical evolution. Let Ni be the predicted total number of long-lived 
stars within the above mentioned subsystem. The predicted counterparts 
for the whole system, related to the selected models, are dA^ = X^i dNi and 
N = J2i Ni, respectively. 

The relative frequency, dN/N, reads: 



dN 



E 



dN, 



E 



A^. dNi 
N N, 



(26) 



A^ 



i=l 



N 
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or, using Eq. I^: 

N M Ni ' ^ ^ 

1=1 ' 

where Mj and M are the mass of the subsystem and the whole system, 
respectively. 

The TDMD related to the whole system, keeping in mind Eq. ( fT5l) . is: 

dN ^ dN, , ^N, dN, 

and the combination of Eqs. ([9]), ( |T5l) . and ( l28i) yields: 

^ = logZ!liTexpio(ai0 + &i) ; (29) 

where the particularization to the thick disk, the thick + thin disk, and the 
Galaxy, is defined by Eqs. f|T3l) and (|T^ . 

Two alternatives shall be considered, Ki and K2 (halo-like + bulge-like 
thick disk), Di and D2 (thick + thin disk), Gi and G2 (Galaxy), according 
if the coefficients, and 6j, are taken from cases HKl, BKl, or HK2, BK2, 
of Table EJ DNl or DN2 of Table El HI, Bl, or H2, B2, of Table 8 (GOT), 
respectively. The last models are related to the halo (H) and bulge (B), 
which are needed in dealing with the Galaxy. 

The TDMD related to the halo-like + bulge-like thick disk and to the thick 
-t- thin disk, is plotted in Fig. [6] (top and bottom panels, respectively) and 
compared to its empirical counterpart, represented in Fig.E] (crosses), both 
in presence (left panels) and in absence (right panels) of [0/Fe] plateau. A 
discontinuity is exhibited by the TDMD, where bulge-like thick disk or thin 
disk formation is assumed to start. 

The TDMD related to the Galaxy is plotted in Fig. [7] and compared to its 
empirical counterpart represented in Fig. H] (middle panels), both in presence 
(left panels) and in absence (right panels) of [0/Fe] plateau. Discontinuities 
are exhibited by the TDMD where bulge, halo-like thich disk, bulge-like thick 
disk, and thin disk formation are assumed to start. 

Although a G-dwarf problem seems to exist for both the halo, the bulge, 
the halo-like thick disk, the bulge-like thick disk, and the thin disk, the 
current model provides a viable interpretation to the occurrence of three 
extremum points, two maxima and one minimum, in the EDMD. 
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5 Inhomogeneous, simple models 



Inhomogeneous simple models with (inhibiting star formation) gas have widely 
been discussed in earlier attempts (COl; C07), following the basic formula- 
tion of the model (Malinie et al., 1993; COO), and the interested reader is 
addressed therein for further details. What is relevant for the current paper, 
shall be reported in the following. 

In the light of the model, each Galactic subsystem is conceived as being 
structured into a number of discrete, entirely gaseous, identical regions, and 
a background of long-lived stars, stellar remnants, and gas inhibited from 
star formation, which have been generated earlier. The evolution occurs via 
a sequence of identical time steps. At the beginning of each step, star forma- 
tion stochastically takes place in a subclass of "active" regions, as described 
by simple homogeneous models, while the others remain "quiescent". At 
the end of each step, high-mass stars have died whereas low-mass stars have 
survived up until today, according to instantaneous recycling approximation. 
In addition, the enriched gas which remains from active regions is instanta- 
neously mixed with the unenriched gas within quiescent regions, to form a 
new set of identical regions for the next step. 

To allow comparison between simple homogeneous and inhomogeneous 
models, the current investigation shall be restricted to the special case of 
expected evolution, where the fraction of active regions is time-independent. 
With respect to the general case, both the theoretical age-metallicity relation 
(TAMR) and the TDMD exhibit minor changes for the disk solar neighbour- 
hood (COO) and negligible variations for the halo solar neighbourhood (COl). 

In any case, active regions evolve as in the Hartwick's (1976) model, 
yielding the usual relationship between the newly synthesised oxygen gas 
mass fraction within an active region, A*(Zq)r, and the gas mass fraction 
within an active region, /j^, at the end of each step in the following way: 

i^*{z'o)K^-p"Hi^'^) ; (30) 

provided that the IMF is universal and the effects of initial composition on 
star evolution are neglected. 

The newly synthesised oxygen gas mass fraction, due to an active region, 
within the system, A*(Zo)r, at the end of each step, is (COO): 

where n is the number of active and quiescent regions, x the probability of a 
region being active, and both remain constant during the expected evolution. 
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The newly synthesised oxygen gas mass fraction, due to all active regions, 
within the system, A*Zo, at the end of each step, is (COO): 



A*Zo = kA*{Zo)R = nxA*{Zo)K ; (32) 

where k = nx is the number of active regions. 

The changes in (allowing star formation) gas and oxygen abundance, re- 
lated to the {i + l)-th step, read (COO; COl): 

fj.i+1 = fieq = fioq'^^^ ; (33) 
(Zo),+i = iZo)e + A*Zo = {Zo)o + (^ + l)A*Zo ; (34) 

g = 1 - X(l - /ik) ; < g < 1 ; (36) 

where the parameter, g, may be considered as an effective gas mass fraction 
within a region at the end of each step, i.e. the mean gas mass fraction 
averaged on both active and quiescent regions. 
The combination of Eqs. flHUl) - flHUl) yields: 



(Zo)£+i = (Zo)o -p'[ln(/X£+i) - ln(/io)] ; (37) 
A*Zo = -p\nq ; (38) 

(1 -/i^)glng 

where p' is an effective yield due to both inhomogeneous star formation and 
inhibited star formation. In the limit of instantaneous recycling i.e. chem- 
ically homogeneous interstellar medium (e.g., the Simple model, Hartwick's 
model), X 1) <? /^r; and p' p" . 

According to Eq. ( 1371) . the dependence of oxygen abundance on the gas 
mass fraction in inhomogeneous models remains unchanged with respect to 
homogeneous models provided that: (i) the expected evolution is considered; 
(ii) the probability of a region being active, xt-i the gas mass fraction within a 
region, /x^^, and the (true) yield, p^, during the (£ + l)-th step, do not depend 
on the evolution; and (iii) the effective yield, p", is replaced by the effective 
yield, p^. In absence of inhibited star formation, it can be shown that the 
effective yield, p', takes the same expression as in the inhomogeneous model 
of Malinie et al. (1993). For further details refer to earlier papers (COO, 
COl). 

Although the TDMD cannot be analytically expressed in the framework 
of inhomogeneous simple models, it still can in the special case of expected 



25 



evolution, with regard to the starting point, ip{4>o) = V'o, and the ending 
point of the initial step, ipi^ipo + ^*(f>o) = 'ipi, where A*0o is the net oxygen 
abundance (normalized to the solar value) increase in gas component at the 
end of the first step. The related explicit expressions can be found in previous 
work [COl; C07, Eqs. (39) and (40) therein]. An approximate expression of 
the TDMD related to the first step, ipi, with the terms up to the second 
order retained, can be found in previous work [COl; C07, Eq. (41) therein]. 
The last expression is valid also in the general case, provided A*(f) is replaced 
by A*(f)o and the probability, by the relative frequency, Uo = ko/rio, being 
ko and Uo the number of active and all regions, respectively, with regard to 
the first step. 

The values of some parameters related to the expected evolution, con- 
cerning cases HKl, HK2; BKl, BK2; DNl, DN2; are listed in Table [71 

The parameters appearing therein, which equal their counterparts corre- 
sponding to homogeneous models (Tables [5] and El cases HK and BK), must 
be connected with active regions, and for this reason some corresponding out- 
put parameters show different values. For further details refer to earlier work 
(COT, Appendix B). It is worth recalling that cases DNl and DN2 listed in 
Table El for comparison, are related to inhomogeneous models listed in Table 

m 

New input parameters are: the effective yield, p', and the normalized 
oxygen abundance increase at the end of a step with regard to the whole 
system, A*0, and to an active region, A*0^, respectively. In presence of 
inhibited (from forming stars) gas (cases HK and BK), the star and inhibited 
gas mass fraction related to the whole system at the end of evolution, Sf and 
Df, are expressed in terms of their counterparts related to an active region 
at the end of a step, and D'^^, as (COT): 

Sf = ^4 ; (40) 

Df = 'j^D'^ ■ (41) 

which are listed in the lower part of Tab. [71 For further details refer to earlier 
work (COT). 

The TDMD related to cases HKl, BKl, DNl, (full lines) and HK2, BK2, 
DN2, (dashed lines) is compared in Fig. [HI to the corresponding EDMD with 
regard to halo- like thick disk (top panels), bulge-like thick disk (middle pan- 
els), and thin disk (bottom panels) stars, both in presence (left panels) and 
in absence (right panels) of [0/Fe] plateau. 

Inhomogeneous simple models with inhibited star formation provide an 
acceptable fit using values of input parameters listed in Table [5l and 0.031 < 
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Table 7: Values of parameters related to the expected evolution of inhomoge- 
neous simple models, in connection with two different cases for the halo-like 
thick disk, HKl and HK2, for the bulge-like thick disk, BKl and BK2, and 
for the thin disk, DNl and DN2, respectively. The indices, 2.9 and 2.35, 
denote values related to the power-law IMF exponent, p, in computing the 
corresponding quantities. For the parameter, ipi, upper and lower values are 
calculated as in C07 by use of Eqs. (40) and (41) therein, respectively. The 
effective yield, p', is related to inhomogenities in oxygen abundance due to 
the presence of active and quiescent regions, whereas oxygen is uniformly 
distributed within active regions. The lower part of the table (last 6 rows) 
is related to models with inhibited star formation. Parameters not reported 
therein have the same value as in the upper part, with the exception of p, 
a, and rrimf, which are listed in Tab. [5] together with other parameters not 
appearing here. The effective yield, p", due to the presence of (inhibiting 
star formation) gas within active regions, is listed as p in the upper part 
of the table. The effective yield, p', due to the presence of both (inhibiting 
star formation) gas within active regions, and (precluding star formation) gas 
within quiescent regions, is listed with the same notation in the upper part of 
the table. The index, R, denotes a generic active region. The mean oxygen 
abundance (normalized to the solar value) of stars at the end of evolution is 
denoted as 0. 
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V' I{Zo)q < 0.039 (halo-like thick disk); 0.094 < p'/{Zo)q < 0.167 (bulge-like 
thick disk); in and p'/{Zo)q ~ 0.442 (thin disk). 

In dealing with active regions, the normalized oxygen abundance at the 
end of evolution, 0/, has to be replaced with the normalized oxygen abun- 
dance at the end of a step, 0j + A*0^. 

Further inspection of Fig.[H] shows that, with respect to homogeneous 
simple models plotted in Fig.El the fit is more or less unchanged for the 
halo-like thick disk, slightly worsened for the bulge-like thick disk in presence 
of [0/Fe] plateau, but substantially improved for the bulge-like thick disk in 
absence of [0/Fe] plateau and the thin disk. In addition, a G-dwarf problem 
still remains (but alleviated) only for the bulge-like thick disk in presence of 
[0/Fe] plateau. 

The TDMD related to a system made of different subsystems, has to 
be numerically computed using Eqs. ( fT3i) and ( l28l) for assigned values of the 
thick disk to thin disk mass ratio, Mkd/^nd, and halo-like to bulge-like 
thick disk mass ratio, Mhk/^bk. The following models have been consid- 
ered with regard to the reference case: Kl and K2 for the thick disk, using 
models HKl, BKl, and HK2, BK2, respectively, defined in Table [71 Dl and 
D2 for the thick + thin disk, using models HKl, BKl, DNl, and HK2, BK2, 
DN2, respectively, defined in Table [71 A comparison with the related EDMD 
(Fig.[3l top and bottom panels, respectively, crosses) is made in Fig.[9l (top 
and bottom panels, respectively) both in presence (left panels) and in ab- 
sence (right panels) of [0/Fe] plateau. The trend looks like its counterpart 
exhibited by homogeneous simple models (Fig. [6]) but the fit is improved, and 
similar considerations can be made. 

The following models have also been considered: Gl and G2 for the 
Galaxy, using models HI, Bl, HKl, BKl, DNl, and H2, B2, HK2, BK2, 
DN2, respectively, where models H and B are related to the halo and the 
bulge, respectively, and are defined in C07 (Table 9). A comparison with the 
related EDMD (Fig.[ll middle and bottom panels, crosses) is made in Fig.[TO] 
both in presence (left panels) and in absence (right panels) of [0/Fe] plateau. 
The trend looks like its counterpart exhibited by homogeneous simple models 
(Fig. [7]) but the fit is improved, and similar considerations can be made. 

In the framework of inhomogeneous models, both active and quiescent 
regins must take origin from (allowing star formation) gas. An inspection of 
Table [71 shows that, for models HK, only a few percent (2.9-5.6%) of the total 
mass has been converted into (long-lived) stars, leaving about 74-39% in (al- 
lowing star formation) gas; and about 22-55% in (inhibiting star formation) 
gas; for models BK, about one quarter (24-27%) of the total mass has been 
converted into (long-lived) stars, leaving about 51-4% in (allowing star for- 
mation) gas; and about 29-69% in (inhibiting star formation) gas; for models 
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DN, more than two thirds (71-87%) of the total mass has been converted into 
(long-hved) stars, leaving about 29-13% in (allowing star formation) gas; and 
no (inhibiting star formation) gas. 

The mean normalized oxygen abundance within long-lived stars is = 

0. 077-0.070, 0.37-0.30, and 0.60-0.50, in the halo-like thick disk, the bulge-like 
thick disk, and the thin disk, respectively. 

It is worth noting that the Galaxy would be mostly luminous for a IMF 
exponent, p = 2.9, and lower stellar mass limit, ^ 0.34, and mostly dark 
for p = 2.35 and mmf ~ 0.007. On the other hand, the assumption of an 
universal IMF implies the same value of the true yield and lower mass limit 
of long-lived stars, for Galactic subsystems. 

6 Discussion 

Although the Galaxy appears to be dominated (by mass and star number) 
by the two major subsystems, the bulge and the thin disk, still the contrary 
holds with respect to selected metallicity bins. For instance, low-metallicity 
bins are dominated by halo stars. Even if the spheroid and the disk evolved 
separately (albeit with similar present metal abundance), the distribution of 
specific angular momentum is known to be consistent with both a halo-bulge 
and a thick disk-thin disk (Wyse and Gilmore, 1992; Ibata and Gilmore, 
1995) collapse. Then in early times a reversed situation occurred, where 
the main Galactic components were the halo and the thick disk, and gas 
exchange between the two reservoirs cannot be excluded. Some signature 
of the formation and the early evolution of the Galaxy could be found in 
the EDMD related to the whole system. In this view, the thick disk has to 
be conceived as pre-existing with respect to the thin disk. For a discussion 
of different scenarios on the origin of the thick disk see e.g., Gilmore et al. 
(1995). 

In brief, the method used in the current and earlier work (C07) is the 
following. An assigned system (e.g., the Galactic spheroid, the Galactic 
thick disk, the Galactic disk, the Galaxy) is conceived as made of a number 
of subsystems where the chemical evolution occurs separately. Then the 
resulting EDMD and TDMD are calculated using Eqs. (E]), ([ED, and (|28l). 

1. e. weighting by mass. The larger uncertainty concerns the thick disk, where 
(i) two different samples (WG95 and CBOO) have been renormalized under 
the assumption that they are equally representative of objects belonging to 
the metallicity range in common, and (ii) the resulting sample (FS07), has 
been conceived as made of two different kind of objects, belonging to the 
halo-like or the bulge-like thick disk, according if the related low-metallicity 
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and the high-metallicity EDMD resembles its halo and bulge counterpart, 
respectively. In this view, the results are merely indicative and need to be 
confirmed using more complete samples. In any case, the G-dwarf problem in 
the bulge-like thick disk, the thick + thin disk, and the Galaxy, is alleviated 
for homogeneous models and much strongly for inhomogeneous models, with 
respect to models where chemical evolution occurs at the same extent in 
different subsystems. 

The absence of a G-dwarf problem in the (thick + thin) Galactic disk 
has been established in a recent attempt (Haywood, 2006), using the Simple 
model. It is not in contradiction with the results of the current paper, for 
the following reasons. First, scalelength corrections in thin disk samples have 
been used in Haywood (2006), while different samples belonging to the thick 
and the thin disk have been used in the current paper. Second, the chemical 
evolution has been assumed to occur at the same rate in the whole disk in 
Haywood (2006), while a different extent has been considered in the current 
paper, according if the halo-like thick disk, the bulge-like thick disk, and the 
thin disk, are dealt with. 

A lower mass limit of long-lived stars, exceeding the theoretical Jeans 
stellar mass (0.007 < nij < 0.01), occurs for a power-law exponent p = 2.9, 
which is a fit to the Scalo (1986) IMF for m ~ m©, concerning both mass 
distribution and oxygen production (Wang and Silk, 1993). A less steep 
Salpeter (1955) IMF, implying p = 2.35, is marginally consistent with the 
theoretical Jeans stellar mass (Tab.E])- On the other hand, the occurrence 
of stellar wind would reduce oxygen nucleosynthesis by a factor of about 2.5 
(e.g., Wang and Silk, 1993) which, in turn, would raise the lower mass limit 
of long-lived stars, provided the (true) yield remains unchanged (e.g., C07, 
Appendix C). 

Within the current model, the lower mass limit of long-lived stars is re- 
lated to a power-law IMF where 2.35 < p < 2.9 provides a poor fit to low 
(m < 1) masses (e.g., Tinney, 1993; Burrows et al., 1993; Binney, 1999; 
Weidner and Kroupa, 2005). In particular, the number of low-mass stars is 
overestimated by a power-law IMF in comparison with current observations. 
The real and the power-law IMF, (f){m) = m~^, can be constrained to yield 
the same mass in a star generation, as: 

/ ^ m<I>reai(m) dm = ; (42a) 

(rnm/)real < ; (42b) 

where rriMf is the upper mass limit of formed stars. Accordingly, the param- 
eter, rrimf, related to a power-law IMF, has to be conceived as an effective 
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lower mass limit of long-lived stars. 

The more relevant parameters of inhomogeneous simple models related to 
halo and bulge (COT), halo-like thick disk, bulge-like thick disk, and thin disk, 
are listed in Table [HI The related values are taken from C07 (Table 9, models 
HI and Bl), and Table [7] (models HKl, BKl, and DNl), respectively. An 
universal power-law IMF is assumed, which makes no change in value of the 
physical parameters, p, a, and m^/. Negative values of inhibition parameter, 
K, and final (inhibiting star formation) gas mass fraction, Df, related to 
the bulge, correspond to enhanced instead of inhibited star formation. For 
further details refer to earlier work (C07). 

For active regions, the probability of star formation is x ~ 1 the 
halo, X < 1 foi' the bulge, and x ^ 1 for both the thick and the thin disk. 
Accordingly, the gas mass fraction remaining at the end of each step, fi'^, is 
close to zero for the halo, about one seventh for the bulge, one fifth-one fourth 
for the thick disk, and one third for the thin disk. The effective yield, p", is 
about twice larger in the thin disk (where it coincides with the true yield) 
than in the bulge-like thick disk; ten times than in the halo-like thick disk; 
five times than in the halo; (where gas is partially inhibited from forming 
stars); and about two thirds lower than in the bulge (where star formation 
is enhanced by inflowing gas with same composition as the preexisting gas; 
for further details refer to C07). 

For the whole system, a similar trend occurs. The ratio of gas mass 
fraction at the end and at the beginning of a step is g <^ 1 for the halo, g < 1 
for the bulge, and g ~ 1 for both the thick and the thin disk. Accordingly, 
the (allowing star formation) gas mass fraction remaining at the end of the 
evolution, /x/, is close to zero for the halo, about four thousandths for the 
bulge, three quarters for the halo-like thick disk, one half for the bulge- 
like thick disk, and one third for the thin disk. The effective yield, p', is 
about three times larger in the thin disk (where it is due to inhomogeneous 
star formation) than in the bulge-like thick disk; a dozen times than in the 
halo-like thick disk; eighty times than in the halo; (where it is due to both 
inhomogeneous star formation and inhibition from forming stars); and about 
three fifths lower than in the bulge (where it is due to both inhomogeneous 
star formation and enhanced star formation; for further details refer to COT). 
At the end of evolution, it remains about 18% of long-lived stars (including 
stellar remnants) and 82% of gas inhibited from forming stars in the halo; 
about 72% of long-lived stars from primeval gas and 28% from inflowed gas 
in the bulge (with respect to the final mass, while mass fractions listed in 
Tableware related to the initial mass); about 0.3% of long-lived stars, 74% 
of (allowing star formation) gas, and 23% of (inhibiting star formation) gas 
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Table 8: Comparison between parameters of inhomogeneous simple models, 
related to the Galactic halo and bulge (CO 7, Table 9, cases HI and Bl), halo- 
like thick disk, bulge-like thick disk, and thin disk (Table [TJ cases HKl, BKl, 
and DNl). An universal power-law initial mass function (IMF) is assumed in 
all cases, which leaves other parameters i.e. p, a, and rrimf, unchanged. Gas 
and star mass fractions, fif, Dj, and Sf, are related to the initial mass with 
regard to the bulge, where total mass conservation is violated by gas inflow 
(COT). Caption of parameters: x - probability of a region being active; /i^ - 
(allowing star formation) gas mass fraction within an active region at the end 
of a step; q - effective gas mass fraction within a region at the end of a step; 
K, - ratio of gas mass fraction which inhibits star formation to long-lived star 
and stellar remnant mass fraction (inhibition parameter); fif - (allowing star 
formation) gas mass fraction at the end of evolution; -D/ - (inhibiting star 
formation) gas mass fraction at the end of evolution; Sf - long-lived star mass 
fraction at the end of evolution; p"/{Zo)q - effective oxygen yield related 
to inhibited or enhanched star formation, normalized to the solar oxygen 
abundance; p'/{Zo)q - effective oxygen yield related to both inhomogeneous 
star formation and inhibited or enhanched star formation, normalized to the 
solar oxygen abundance. 



parameter value 





halo 




bulge 




halo-like 
thick disk 




bulge-like 
thick disk 




thin disk 




X 


9.5579 10" 


-1 


8.7545 10" 


-1 


9.3147 10" 


-2 


1.9961 10- 


-1 


7.6765 10- 


-2 




8.2201 10" 


-4 


1.4139 10" 


-1 


2.3356 10- 


-1 


2.1579 10- 


-1 


3.7148 10- 


-1 


Q 


4.4991 10- 


-2 


2.4833 10- 


-1 


9.2861 10- 


-1 


8.4346 10- 


-1 


9.5175 10- 


-1 


K 


4.5169 




-3.2099 10- 


-1 


7.8271 




1.0370 




0.0000 




fif 


4.0973 10" 


-6 


3.8030 10" 


-3 


7.4358 10- 


-1 


5.0613 10- 


-1 


2.9047 10- 


-1 


Df 


8.1874 10" 


-1 


-4.7094 10" 


-1 


2.2737 10- 


-1 


2.5142 10- 


-1 


0.0000 






1.8126 10" 


-1 


1.4671 




2.9049 10- 


-2 


2.4245 10- 


-1 


7.0953 10- 


-1 




1.3363 10" 


-1 


1.0857 




8.3518 10- 


-2 


3.6191 10- 


-1 


7.3722 10- 


-1 


P'/iZo)e 


5.3452 10~ 


-3 


7.6001 10- 


-1 


3.8418 10- 


-2 


1.6648 10- 


-1 


4.4233 10- 


-1 
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in the halo-like thick disk; about 24% of long-lived stars, 51% of (allowing 
star formation) gas, and 25% of (inhibiting star formation) gas in the bulge- 
like thick disk; about 71% of long-lived stars and 29% of gas allowing star 
formation in the thin disk. 

The TDMD related to both homogeneous and inhomogeneous simple 
models provides an acceptable fit to the EDMD related to the thick disk, 
the thick + thin disk, (Figs.E] and upper and lower panels, respectively) 
and the Galaxy (Figs. [7] and [TOl) : in particular, a non monotonic trend is 
reproduced. While models assume that star formation at the beginning of 
the evolution of each subsystem starts abruptly with constant efficiency, the 
data seem to indicate a somewhat gradual rate, with increasing efficiency. 
In other words, the occurrence of some physical process (not necessarily gas 
infall) seems to inhibit star formation in early times of evolution. Accord- 
ingly, the history of each Galactic subsystem (halo, bulge, halo-like thick disk, 
bulge-like thick disk, thin disk) could be conceived as made of two distinct 
phases, namely (i) assembling, where (long-lived) star formation efficiency is 
gradually increasing, and (ii) stabilization, where (long-lived) star formation 
efficiency maintains constant. 

The assumption that the EDMD of the local disk is representative of 
the global disk, even if in contrast with an inside-out disk formation, can 
be considered as a useful zero-th order approximation. On the other hand, 
nearby stars older than about 0.2 Gyr come from birth sites which span a 
large range in Galactocentric distances (e.g., Rocha-Pinto et al., 2000). The 
orbital diffusion coefficient deduced from the observed increase of velocity 
dispersion with age implies that presently local stars have suffered a rms 
azimuthal drift from about 2 kpc (for an age of 0.2 Gyr) to many Galactic 
orbits (for an age of 10 Gyr); for further details refer to earlier work (Wie- 
len, 1977). Considerable, but smaller, drift should occur also on the radial 
direction. In this sense, the star formation rate inferred for nearby stars is 
a measure of the global Milky Way star formation rate, at least at the sun 
Galactocentric radius (Rocha-Pinto et al., 2000), according to the estimates 
of the diffusion coefficient (e.g., Meusinger et al., 1991). 

As outlined in an earlier attempt (Wang and Silk, 1993), radial flows of 
the gas (possibly due to the observed spiral density waves) or different disk 
star formation histories (with time scales comparable to that of chemical 
evolution) between the inner end the outer parts of the Galaxy, may change 
the local abundances, but the overall abundance in the disk should not be 
affected as long as there is no gain or loss of material in the disk. It can 
also be noticed (Wang and Silk, 1993) that the average oxygen abundance 
in the disk is roughly the solar value, and can be plausibly explained by the 
standard Scalo IMF with lower and upper star mass limit, m^/ = 0.1 and 



33 



TriMf = 60, respectively, or by a power-law IMF with exponent, p, within the 
range, 2.35 < p < 2.9, and the same rfimf and niMf- For further details refer 
to the parent paper (Wang and Silk, 1993). 

At the end of halo-like and bulge-like thick disk evolution, the fractional 
gas and star mass predicted by the model are: 



. = (/^uk)/ + {DvK)f ; U = H, B ; (43) 

U^^UK Jo 



(M, 

{M 
(Muk) stars 



UKjo 



(suk)/ ; U = H,B ; (44) 



where (Mhk)o and (Mbk)o are the initial halo-like and bulge-like thick disk 
mass, respectively. The combination of Eqs. 0431) and (jH]) yields: 

(MuK)gas = ^^^^4^^^^^(MuK)stars U = H, B ; (45) 

where (MHK)stars and (MBK)stars must be equal to the current halo-like and 
bulge-like thick disk mass, respectively. 

For models HKl-2, BKl-2, the values of fractional masses, (/xuk) /, (-Duk) /, 
and (suk)/, are listed in Table [71 and taking a current halo- like thick disk 
mass, Mhk = (MHK)stars = (9/296) (29/55)Mio ^ O.OI6M10, and a current 
bulge-like thick disk mass, Mbk = (MBK)stars = (287/296) (29/55)Mio ^ 
O.5IIM10, according to Eqs. fll3bp and ffT^ . it is found the following: 



(MHK)gas = 0.5359 - 0.2705Mio 
(MBK)gas = 1.5974 - I.335IM10 
(MDK)gas = 2.1333 - 1.6056Mio 



(46) 
(47) 
(48) 



where (MDK)gas = (^HK)gas + (^BK)gas is the thick disk gas mass budget 
available for the thin disk formation, which amounts to about 36%-28% of 
the assumed thin disk mass, Eq. ( ]13b[) . 

The above discrepancy could be avoided, taking into consideration mass 
loss during the thick disk evolution, in the sense that gas flows away but 
chemical evolution therein occurs as in the remaining gas. Accordingly, 
Eqs. (H3l)- (H5l) still hold but the star mass at the end of halo-like and bulge- 
like thick disk evolution is related to current halo-like and bulge-like thick 
disk mass, as: 

(MuK)stars = ; U = H, B ; (49) 

1 - Ju 

where (MuK)stars includes stars both inside and outside the thick disk, and 
/h, /b, are the star mass fraction outside the halo-like and bulge-like thick 
disk, respectively. 
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The boundary conditions to the problem under discussion are: 



(MHK)gas = Mbk + Mnd ; (50) 

(MBK)gas = MND ; (51) 

the gas mass budget at the end of the halo-hke and bulge-hke thick disk must 
necessarily equal the bulge-like thick disk + thin disk and the thin disk mass, 
respectively, in absence of mass loss outside the disk. 

The combination of Eqs. (US]), (HHD, dnUD, and ([511) yields: 

^huMbk + M^,D _ ^ (/.uk)/ + {DvK)f . y = H, B ; (52) 

where 6ij is the Kronecker symbol. Then the halo-like and bulge-like thick 
disk mass fraction to be lost during the evolution is: 

^ ^ _ ijMjK}i+iD^ ^ ; U = H,B ; (53) 

(suk) / OuvMbk + Mnd 



which, for models HKl-2, BKl-2, and mass values prescribed by Eqs. (113bp 
and f|T^ . reads: 

/jj = 0.9074 - 0.9532 ; (54) 
/b = 0.6970 - 0.7468 ; (55) 

where it can be seen that more than 90% of the halo-like thick disk initial 
gas budget and more than 69% of the bulge-like thick disk initial gas budget 
were removed from their reservoirs during the evolution. In other words, the 
bulge-like thick disk and the thin disk did not begin to start abruptly, but 
were gradually assembled during the evolution of their precursors, as inferred 
from the features of the EDMD mentioned above. 

The results found in the current paper rely on the assumption of universal 
IMF, but the method used here can be generalized to the case where different 
IMFs in space and/or time are considered. To this aim, it is sufficient to 
realize that the relation, Mj = TrTiNi, with regard to i-th subsystem, holds 
for different mean stellar masses, mj, related to different IMFs. Accordingly, 
Eq. ([9]) has to be generalized as Ni/N = {rn/mi){Mi/ M), where m is the 
mean stellar mass within the system. Different IMFs during the assembling 
and stabilization phase have been considered since a long time (e.g., Schmidt, 
1963; Caimmi, 1978a,b, 1981, 1982; COO). 

A different IMF in the bulge, with respect to the disk, has recently been 
advocated for explaining the giant star metal distribution in both the Milky 
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Way and Andromeda bulge (Ballero et al., 2007a), contrary to what has been 
found for the Milky Way bulge (C07). This discrepancy could be owing to 
two orders of reasons, namely: (i) instantaneous mixing has been assumed 
for chemical evolution in the former case and instantaneous recycling in the 
latter and, (ii) the cumulative metallicity distribution has been used in the 
former case and the differential metallicity distribution in the latter. In ad- 
dition, different imput parameters have been used for the related models. 
A comparison between observations and model predictions with regard to 
[a/Fe] ratios (e.g., Ballero et al., 2007b), cannot be performed in the frame- 
work of the model used here, as only oxygen chemical evolution is considered. 
Neverthless, it is worth mentioning that Fe production mainly depends on 
type la supernovae which, in turn, have binary star progenitors. Then dif- 
ferent IMFs in different subsystems could be mimiced by different amounts 
of binary star fractions. 

The results found in the current paper also rely on the assumed [0/H]- 
[Fe/H] relation, Eqs. ([I])-([2]). Data from recent investigations, different for 
different Galactic subsystems (e.g., Bensby et al., 2004; Jonsell et al., 2005; 
Fulbright et al., 2005; Garcia Perez et al., 2006; Melendez et al., 2006; 
Ramirez et al., 2007) yield [0/H]-[Fe/H] relations lying between their coun- 
terparts expressed by Eqs. ([T]) and ([2]) which, in turn, allow comparison with 
earlier results related to the halo (COl), the bulge (C07), and the Galactic 
spheroid (C07). On the other hand, the method used in the current paper 
holds regardless from the assumed (or not) [0/H]-[Fe/H] relation. 

Although direct [0/H] values are available from recent observations, still 
they cannot be used for the global EDMD in absence of their counterparts 
related to each Galactic subsystem. For instance, a thick disk sample inves- 
tigated by Reddy et al. (2006) is biased towards low-metallicity objects sim- 
ilarly to the WG95 sample; a sample investigated by Ramirez et al. (2007), 
hereafter quoted as the RA07 sample, is made of = 523 nearby stars be- 
longing to the halo, the thick disk, and the thin disk, within a metallicity 
range, —1.5 < [Fe/H]< 0.5. Distinct and homogeneous samples representa- 
tive of the halo, the bulge, the thick disk, and the thin disk, where oxygen 
abundance is directly inferred, would considerably improve the results. 

Neverthless, an indicative comparison between the EDMD deduced for 
the thick -|- thin disk can be made, with regard to the FS07 sample (Fig.[3l 
bottom panels, crosses), and the RA07 sample where oxygen abundance was 
directly inferred both in presence and in absence of the LTE approximation 
(Ramirez et al., 2007, Table 6, available at the CDS). The local thermody- 
namic equilibrium (LTE) approximation well holds deep in the stellar pho- 
tosphere, and a continuous transition occurs up to complete non-equilibrium 
(NLTE) high in the atmosphere. For further details refer to specialistic 
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Table 9: The empirical, differential metallicity distribution (EDMD), de- 
duced from the RA07 sample {N = 523), both in presence (LTE) and in 
absence (NLTE) of the LTE approximation. 
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textbooks (e.g.. Gray, 2005, Chap. 6). The related EDMD, obtained by 
use of Eqs. (l3])-(l6l), is listed in Table M both in presence (LTE) and in ab- 
sence (NLTE) of the LTE approximation. Halo contamination may safely be 
neglected except at low metallicities, where the global contribution is also 
negligible in that (i) the halo star subsample in the RA07 sample is by far 
incomplete, and (ii) the metal-weak thick disk ([Fe/H]< —1) has not been 
considered therein. For further details refer to the parent paper (Ramirez et 
al., 2007). 

The EDMD deduced from the RA07 sample and listed in Table [9l is 
compared in Fig.[TT] with its counterpart related to the thick -|- thin disk, 
plotted in Fig. El (bottom panels, crosses), both in presence and in absence of 
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the LTE approximation in the former case, and [0/Fe] plateau in the latter. 
The low-metallicity tail is zoomed in bottom panels. It is apparent that 
using the LTE approximation understimates the EDMD for subsolar oxygen 
abundances and overstimates for supersolar oxygen abundances. 

Although the EDMD related to the thick + thin disk (Fig. El bottom 
panels, crosses) is systematically lower for low oxygen abundances and higher 
for intermediate values (ending at ~ 1.56) with respect to the EDMD 
deduced from the RA07 sample both in presence (triangles) and in absence 
(diamonds) of the LTE approximation, still a similar trend is exhibited. 

A closer agreement could occur in the low-metallicity range after removing 
halo contamination in the RA07 sample or, alternatively, what in the current 
paper has been defined as the "halo-like thick disk" could be due to halo 
contamination in the CBOO sample. In any case, the RA07 sample cannot be 
considered as representative of the metal-weak ([Fe/H]~ —1) thick disk, as 
mentioned above. The discrepancy in the intermediate-metallicity range is 
reduced in absence of the LTE approximation, and further (shghtly) reduced 
in presence of [Fe/H] plateau (left panels), with the exception of the more 
metal-rich (0 ^ 1.56) bin in the disk EDMD, which is related to the more 
metal-rich bin in the RM96 sample. Accordingly, the RM96 sample could 
be biased towards high metallicities, and the above mentioned bin have little 
statistical meaning (e.g., Haywood, 2001). In conclusion, the disk EDMD 
plotted in Fig.O (bottom panels) may be considered as acceptable to a first 
extent. 

The empirical age-metallicity relation (EAMR), taken from a subsample 
(A^ = 223) of the RA07 sample, is plotted in Figs.fT2] and fT3] in presence and 
in absence of the LTE approximation, respectively. Also represented therein 
is the temporal behaviour of (allowing star formation) gas oxygen abundance, 
normalized to the solar value, related to models HI (C07, Table 9; dashed 
lines), HKl (Tabled dotted lines), BKl (dash-dotted lines), and DNl (full 
lines). The curves correspond to the special case of constant star formation 
rate within active regions but, in any case, the extreme points related to the 
beginning and the end of a step, are left unchanged. Then the curves may 
be considered as representative of a time-dependent star formation rate. 

It is worth recalling that simple models of chemical evolution used in the 
current paper are time- independent, with regard to the TDMD. Accordingly, 
the initial and the final time, together with the time step in the case of 
inhomogeneous models, can be selected for best fitting the EAMR. 

For the halo, according to C07, initial and final values, ([0/H], T/Gyr) = 
(—3, 12.5) and (0, 8.0), respectively, have been chosen, together with a time 
step, AT/Gyr= 1.125, which implies four steps during the evolution. The ini- 
tial and final values of fractional oxygen abundance, {(pe, A*0^), at any step. 
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are (0.001,0.9502), (0.0176,0.9668), (0.0342,0.9834), (0.0508,1), correspond- 
ing to constant increments, A*(j) = (pi — (pe-i = 0.01658 and A*(j)'^ = 0.9492. 
The related TAMR is represented by four (the first out of scale) dashed lines 
in Figs. [12] and [131 Owing to the logarithmic scale on the horisontal axis, 
log0 = [0/H], the ending point of each step is very close to [0/H] = 0, but 
only the last step actually ends therein. The starting point of each step is 
out of scale on the left. 

For the halo-like thick disk, initial and final values, ([0/H], T/Gyr) = 
(—1.52, 11.0) and (—0.80, 10.2), respectively, have been chosen, together with 
a time step, AT/Gyr= 0.2, which implies four steps during the evolution. 
The initial and final values of fractional oxygen abundance, (0^, A*(j)'^), at any 
step, are (0.03,0.1515), (0.03285,0.1543), (0.03569,0.1572), (0.03854,0.16), 
corresponding to constant increments, A*0 = (pi — (pe-i = 0.002846 and 
A*0^ = 0.1215. The related TAMR is represented by four dotted lines in 
Figs.[T2]and[T3l The starting point of each step is out of scale on the left. 

For the bulge-like thick disk, initial and final values, ([0/H], T/Gyr) = 
(—0.80, 10.2) and (—0.61, 9.0), respectively, have been chosen, together with 
a time step, AT/Gyr= 0.3, which implies four steps during the evolution. 
The initial and final values of fractional oxygen abundance, (0£,A*0^), at 
any step, are (0.16,0.7150), (0.1883,0.7433), (0.2167,0.7717), (0.2450,0.8), 
corresponding to constant increments, A*0 = (pi — (pe-i = 0.002834 and 
A*(/)^ = 0.5550. The related TAMR is represented by four dot-dashed lines 
in Figs. [12] and [121 The starting point of the first step coincides with the 
ending point of the last step related to the halo-like thick disk. 

For the thin disk, initial and final values, ([0/H], T/Gyr)= (-0.49, 10.0) 
and (0.20, 0.0), respectively, have been chosen, together with a time step, 
AT/Gyr= 0.4, which implies twenty-five steps during the evolution. The 
initial and final values of fractional oxygen abundance, {(pi, A*0^), at the first 
five steps, are (0.325, 1.0550), (0.3469, 1.0769), (0.3687, 1.0988), (0.3906, 1.1207), 
(0.4125, 1.1425), corresponding to constant increments, A*(p = (pi — (pi-i = 
0.02187 and A*(p'^ = 0.7300. The related TAMR is represented by five full 
lines in Figs. [12] and [131 The remaining twenty steps are out of scale on the 
bottom. 

It is apparent that the TAMR related to the disk cannot provide an 
acceptable fit to the EAMR deduced from the RA07 sample both in presence 
and in absence of the LTE approximation. This is why models of chemical 
evolution used in the current paper have been selected to fit the FS07 sample 
for the thick disk and the RM96 sample for the thin disk, where a small 
amount of oxygen-rich (supersolar) objects were found using Eqs. ([1]) and 
([2]). On the contrary, a consistent number of oxygen- rich objects is present 
in the RA07 sample which, on the other hand, is biased towards oxygen-poor 



39 



objects (Ramirez et al., 2007). 

In the future, additional effort should be devoted in collecting represen- 
tative samples of both the halo, the bulge, the halo-like and bulge-like thick 
disk, and the thin disk, where the oxygen abundance is directly derived (e.g., 
Ramirez et al., 2007) instead of being empirically deduced from iron abun- 
dance, as in Eqs. ([T]) and ([2]). Additional support to the distinction between 
halo-like and bulge-like thick disk could be provided from the related dis- 
tributions of specific angular momentum. In particular, the former and the 
latter are expected to exhibit differences as in the Galactic spheroid with 
respect to the Galactic disk (Wyse and Gilmore, 1992; Ibata and Gilmore, 
1995). It would also be important to establish if the metal- weak thick disk 
(e.g., Martin and Morrison, 1998; Beers et al., 2002) is the actual precur- 
sor of the metal-enriched thick disk or, on the other hand, coincides with 
the flattened, rotation-supported halo subsystem (Prochaska et al., 2000), or 
both. Further investigation will provide additional support, or contradict, 
the results of the current attempt. 

7 Conclusion 

The empirical differential metallicity distribution (EDMD) in the thick disk 
has been deduced from two different samples involving (WG95) local thick 
disk stars derived from Gliese and scaled in situ samples within the range, 
-1.20 <[Fe/H]< -0.20 (Wyse and Gilmore, 1995), and (CBOO) = 46 
likely metal-weak thick disk stars within the range, —2.20 <[Fe/H]< —1.00 
(Chiba and Beers, 2000). The EDMD in the thin disk has been deduced 
from (RM96) = 287 chemically selected G dwarfs within 25 kpc from the 
sun, with the corrections performed in order to take into account the stellar 
scale height (Rocha-Pinto and Maciel, 1996). To this aim, two alternative 
[0/H]-[Fe/H] dependences have been used, Eqs. ([I]) and ([2]), according to 
earlier attempts where the EDMD was deduced in halo field stars (COl), the 
globular cluster subsystem (C07), and the bulge (C07). 

The metal-poor and metal-rich EDMD related to the thick disk, have 
shown analogies with their halo and bulge counterparts, respectively. For this 
reason, the thick disk has been conceived as made of two distinct regions: the 
halo-like and bulge-like thick disk, and the related EDMD has been deduced. 
The data have been fitted, to an acceptable extent, by both homogeneous and 
inhomogeneous simple models of chemical evolution. For the thin disk, only 
inhomogeneous models have suceeded in reproducing the data, as already 
shown in earlier attempts (Malinie at al., 1993; COO). 

Under the assumption of an universal initial mass function (IMF) and a 
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value of the true yield equal to the linear fit related to the thin disk solar 
neighbourhood, inhibition of thick disk star formation (implying gas outflow) 
has been required to reproduce the EDMD. A power-law IMF has been con- 
sidered, 0(m) oc m~^, within the range, 2.35 < p < 2.9, and special effort 
has been given to the limiting cases, (i) p = 2.9, which is acceptably close 
to Scalo (1986) IMF for m ~ m©, and (ii) p = 2.35, which is the Salpeter 
(1955) IMF. 

The EDMD related to the disk has been determined by weighting the 
mass, where each distribution has been assumed to be typical for the corre- 
sponding subsystem. Following a similar procedure, the EDMD related to 
the Galaxy has been computed using the results of an earlier attempt (C07) 
with regard to the bulge and the halo. In any case, it has been inferred 
that a more refined model involving an initially increasing star formation 
efficiency (but not necessarily implying gas infall) while assembling Galactic 
subsystems, could provide a better agreement with the data, including the 
right amount of mass budget needed for the formation of the bulge-like thick 
disk and the thin disk. The theoretical differential metallicity distribution 
(TDMD) related to both homogeneous and inhomogeneous simple models 
has provided an acceptable fit to the EDMD related to the Galaxy and, in 
particular, a non-monotonic trend has been reproduced. 

An indicative comparison has been performed between the EDMD de- 
duced for the thick + thin disk both in presence and in absence of [0/Fe] 
plateau, and its counterpart computed for (RA07) = 523 nearby {d <150 
pc) stars with metallicity in the range, —1.5 <[Fe/H] < 0.5, for which the 
oxygen abundance has been determined both in presence and in absence 
of the LTE approximation (Ramirez et al., 2007). It has been found that 
both distributions exhibit a similar trend, though systematic differences ex- 
ist, which has made the fit acceptable to a first extent. 

The empirical age-metallicity distribution (EAMR), taken from a subsam- 
ple (A^ = 223) of the RA07 sample for which age determination is available 
(Ramirez et al., 2007), has been compared with the theoretical age-metallicity 
relation (TAMR), predicted by models related to the halo, the halo-like thick 
disk, the bulge-like thick disk, and the thin disk. The fit has shown to be 
unsatisfactory for the following reason. The models of chemical evolution 
used in the current paper have been selected to match the data from earlier 
samples for both the thick disk (WG95 and CBOO) and the thin disk (RM96), 
where a small amount of oxygen-rich (supersolar) objects were found using 
Eqs. ([1]) and ([2]). Conversely, a consistent number of oxygen- rich objects has 
been found in the RA07 sample which, on the other hand, is biased towards 
oxygen-poor objects (Ramirez et al., 2007). 
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Figure 1: The empirical, differential metallicity distribution (EDMD) in thick; 
(top panels) and thin (bottom panels) disk, both in presence (PP) and in 
absence (AP) of [0/Fe] plateau. Plots have been deduced from the FS07 
sample for the thick disk, and from the RM96 sample for the thin disk. The 
dotted vertical line marks the transition from halo to bulge/disk globular 
cluster morphological type, [Fc/H]= —0.8. The dashed vertical line marks 
the transition from halo-like to bulge-like thich disk, [Fe/H]= —1.4. 
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Figure 2: The empirical, differential metallicity distribution (EDMD) in halo- 
like (top panels) and bulge-like (bottom panels) thick disk, both in presence 
(PP) and in absence (AP) of [0/Fc] plateau. The dotted vertical line marks 
the transition from halo to bulge/disk globular cluster morphological type, 
[Fe/H]= -0.8. 4g 
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Figure 3: The empirical, differential metallicity distribution (EDMD) in halo- 
like + bulge-hke thick disk (top panels) and halo-like -|- bulge-hke thick -|- 
thin disk (bottom panels), both in presence (PP) and in absence (AP) of 
[0/Fe] plateau. Different symbols are related to different halo-like to bulge- 
like thick disk mass ratio, Mhk/^bk, for top panels, and to different thick 
to thin disk mass ratio, Mkd/A^nd, for bottom panels. Caption of symbols: 
top panels, Mhk/-/Wbk = 18/574 (crosses); 0.01 (triangles); 0.10 (squares); 
bottom panels, Mkd/A^nd = 0.10 (crosses); 0.05 (triangles); 0.20 (squares). 
In addition, Mrk/Mbk = 18/574 ^ 0.03136 for bottom panels. The dotted 
vertical line marks the transition from halo to bulge/disk globular cluster 
morphological type, [Fe/H]= —0.8. 
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Figure 4: The empirical, differential metallicity distribution (EDMD) in 
the Galaxy (middle panels) for halo-like to bulge-hke thick disk mass ratio, 
A^hk/^bk = 18/574, and thick disk to thin disk mass ratio, Mkd/^nd = 
0.1, both in presence (PP) and in absence (AP) of [0/Fe] plateau. The 
EDMD in the Galactic spheroid (C07) is shown in top panels for compar- 
ison. Additional cases related to different mass ratios are also represented 
(bottom panels) where the error bars have been suppressed to avoid con- 
fusion. Caption of symbols: Mhk/^bk = 0.01 (triangles), 0.10 (squares), 
where A^kd/Mnd = 0.1; Mkd/Mnd = 0.05 (diamonds), 0.20 (asterisks), 
where Mhk/-^bk = 18/574. Crosses representing error bars arc related to 
the reference case, Mhk/Mbk = 18/574 ^ 0.03136 and Mkd/Mnd = 0.10. 
The dotted vertical hue marks the transition from halo to bulge/disk globular 
cluster morphological type, [Fe/H]= —0.8. 
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Figure 5: Comparison between theoretical (TDMD) and empirical differential 
metallicity distribution (EDMD) in halo-like thick disk (top panels), bulge- 
like thick disk (middle panels) and thin disk (bottom panels) stars, both 
in presence (left panels) and in absence (right panels) of [0/Fe] plateau, 
respectively. The straight lines correspond to homogeneous models, HKl- 
HK2 (top panels), BK1-BK2 (middle panels), and to inhomogeneous models, 
DN1-DN2 (bottom panels), defined in Table El The dashed line in top panels 
(model HK2) is very close to the limiting case, fif = 0, defined by Eqs. ( |T8l) 
and (fT9l) and model HKl is represented by the less inclined (in absolute 
value) full line. Accordingly, the linear fits lying between the dashed line and 
the more inclined (in absolute value) full hne, must be excluded. Crosses 
represent the data and related uncertainties, as in Figs.[T]|2l The dotted 
vertical line marks the transition from halo to bulge/disk globular cluster 
morphological type, [Fe/H]=— 0.8. 
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Figure 6: Comparison between theoretical (TDMD) and empirical differ- 
ential metallicity distribution (EDMD) in halo-like + bulge-like thick disk 
(top panels) and thick -|- thin disk (bottom panels), both in presence (left 
panels) and in absence (right panels) of [0/Fe] plateau, respectively. The 
straight lines correspond to homogeneous models, HK1-HK2 (low oxygen 
abundances), BK1-BK2 (intermediate oxygen abundances), and DN1-DN2 
(high oxygen abundances) , defined in Table El Crosses represent renormal- 
ized data and related uncertainties, as in Fig.El The dotted vertical line 
marks the transition from halo to bulge/disk globular cluster morphological 
type, [Fe/H]=— 0.8. The dashed vertical line marks the (assumed) transition 
from halo-like to bulge-like thick disk. 
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Figure 7: Comparison between theoretical (TDMD) and empirical differential 
metallicity distribution (EDMD) in the Galaxy (top panels) and zoomed for 
low oxygen abundance (bottom panels), both in presence (left panels) and in 
absence (right panels) of [0/Fe] plateau, respectively. The curves correspond 
to models HI, Bl, HKl, BKl, DNl (full); H2, B2, HK2, BK2, DN2 (dashed) 
defined in: COT, Tables 7 and 8 therein (H, B); Tables [5] and [6] (HK, BK, 
DN); all combined via Eq. (l29ll . The dotted vertical line marks the transition 
from halo to bulge/disk globular cluster morphological type, [Fe/H]=— 0.8. 
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Figure 8: Comparison between theoretical (TDMD) and empirical differ- 
ential metallicity distribution (EDMD) in halo-like thick disk (top panels), 
bulge-like thick disk (middle panels), and thin disk (bottom panels) stars, 
both in presence (left panels) and in absence (right panels) of [0/Fe] plateau, 
respectively. Full curves correspond to inhomogeneous models HKl (top pan- 
els), BKl (middle panels), and DNl (left bottom panel), and dashed curves 
to HK2 (top panels), BK2 (middle panels), and DN2 (right bottom panel), 
defined in Table [71 Crosses represent the data and related uncertainties, 
as in Figs. [mil The dotted vertical line marks the transition from halo to 
bulge/disk morphological type in globular clusters, [Fe/H]=— 0.8. 
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Figure 9: Comparison between theoretical (TDMD) and empirical differen- 
tial metallicity distribution (EDMD) in halo-like + bulge-like thick disk (top 
panels) and thick -|- thin disk (bottom panels), both in presence (left panels) 
and in absence (right panels) of [0/Fe] plateau, respectively. Full and dashed 
curves correspond to models Kl, K2, (top panels), and Dl, D2, (bottom pan- 
els), obtained using models HKl, BKl; HK2, BK2; (top panels); and HKl, 
BKl, DNl; HK2, BK2, DN2; (bottom panels); respectively, defined in Tables 
[5] and [71 and combined via Eqs. ([9]) and (!28|) . Crosses represent the renor- 
malized data and related uncertainties, as in Fig. [HI The dotted vertical line 
marks the transition from halo to bulge/disk globular cluster morphological 
type, [Fe/H]=— 0.8. The dashed vertical line marks the (assumed) transition 
from halo-like to bulge-like thick disk. 
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Figure 10: Comparison between theoretical (TDMD) and empirical differen- 
tial metallicity distribution (EDMD) in the Galaxy (top panels) and zoomed 
for low normalized oxygen abundances (bottom panels), both in presence 
(left panels) and in absence (right panels) of [0/Fe] plateau, respectively. 
Curves correspond to models HI, Bl, HKl, BK1,DN1, (full); H2, B2, HK2, 
BK2, DN2, (dashed); defined in C07, Tables 7 and 9 therein (H, B) and Ta- 
bles [5] and [7] (HK, BK, DN), all combined via Eqs. (El) and (EHl). The dotted 
vertical line marks the transition from halo to bulge/disk globular cluster 
morphological type, [Fe/H]=— 0.8. 
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Figure 11: The empirical, differential metallicity distribution (EDMD) de- 
duced from the RA07 sample (A^ = 523) including nearby halo (by far in- 
complete), thick disk, and thin disk stars, both in presence (triangles) and in 
absence (diamonds) of the LTE approximation. The EDMD deduced for the 
thick -|- thin disk (Fig.[3l crosses) is also represented for comparison both in 
presence (left panels) and in absence (right panels) of [0/Fe] plateau. The 
low-met allicity tail is zoomed in bottom panels. 
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Figure 12: Comparison between the empirical (EAMR) and theoretical age- 
metallicity relation (TAMR) in presence of the LTE approximation. The data 
are from a subsample {N = 223) of the RA07 sample, for which ages have 
been determined (Ramirez et al., 2007). Dashed, dotted, dash-dotted, and 
full lines arc related to models HI (C07), HKl, BKl, and DNl, respectively, 
in the special case of constant star formation rate within active regions. 
The halo star formation begins at ([0/H], T/Gyr)= (—3, 12.5) and ends at 
(0, 8.0), within four time steps, AT/Gyr=1.125 (C07). The starting point of 
each step is out of scale on the left; the first step is out of scale on the top; and 
cannot be shown. The halo-like thick disk star formation begins at (—1.52, 
11.0) and ends at (-0.80, 10.2), within four time steps, AT/Gyr=0.2. The 
starting point of each step is out of scale on the left, and cannot be shown. 
The bulge-like thick disk star formation begins at (—0.80, 10.2) and ends 
at (—0.10, 9.0), within four time steps, AT/Gyr=0.3. The thin disk star 
formation begins at (—0.49, 10.0) and is still continuing (in the sense that 
star formation occurs) up today at (0.20, 0.0) within twenty-five time steps, 
AT/ Gyr=0.4. The last twenty time steps are out scale on the bottom. 
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Figure 13: Comparison between the empirical (EAMR) and theoretical age- 
metallicity relation (TAMR) in absence of the LTE approximation. The data 
are from a subsample {N = 223) of the RA07 sample, for which ages have 
been determined (Ramirez et al., 2007). Dashed, dotted, dash-dotted, and 
full lines arc related to models HI (C07), HKl, BKl, and DNl, respectively, 
in the special case of constant star formation rate within active regions. 
The halo star formation begins at ([0/H], T/Gyr)= (—3, 12.5) and ends at 
(0, 8.0), within four time steps, AT/Gyr=1.125 (C07). The starting point of 
each step is out of scale on the left; the first step is out of scale on the top; and 
cannot be shown. The halo-like thick disk star formation begins at (—1.52, 
11.0) and ends at (-0.80, 10.2), within four time steps, AT/Gyr=0.2. The 
starting point of each step is out of scale on the left, and cannot be shown. 
The bulge-like thick disk star formation begins at (—0.80, 10.2) and ends 
at (—0.10, 9.0), within four time steps, AT/Gyr=0.3. The thin disk star 
formation begins at (—0.49, 10.0) and is still continuing (in the sense that 
star formation occurs) up today at (0.20, 0.0) within twenty-five time steps, 
AT/ Gyr=0.4. The last twenty time steps are out scale on the bottom. 
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